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Abstract: In these lectures I give a pedagogical presentation of some of the recent progress in 
supersymmetric Chern-Simons-matter theories, coming from the use of localization and matrix 
model techniques. The goal is to provide a simple derivation of the exact interpolating function 
for the free energy of ABJM theory on the three-sphere, which implies in particular the iV 3 / 2 
behavior at strong coupling. I explain in detail part of the background needed to understand this 
derivation, like holographic renormalization, localization of path integrals, and large ./V techniques 
in matrix models. 
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1. Introduction 

The AdS / CFT correspondence postulates a remarkable equivalence between certain gauge theo- 
ries and string/M-theory on backgrounds involving AdS spaces. In these lectures we will look at 
the correspondence in the case of AdS4/CFT3, focusing on ABJM theory [6]. One of the myste- 
rious consequences of this correspondence is that, at strong coupling, the number of degrees of 
freedom of the CFT3 should scale as TV 3 / 2 [68] . This was finally established at the gauge theory 
level in [32]. 

The derivation of the TV 3 / 2 behavior in [32] is based on various ingredients. The first one, due 
to [61], is the fact that certain path integrals in superconformal Chern-Simons-matter theories 
can be reduced to matrix integrals, by using the localization techniques of [82]. Localization 
techniques have a long story in supersymmetric and topological QFTs, and the applications of [82, 
61] to superconformal field theories provide a powerful technique to analyze certain observables 
in terms of matrix models. One of these observables is, in the three-dimensional case, the free 
energy on the three-sphere. In [32] it was pointed out that this quantity is a good measure of the 
number of degrees of freedom in a CFT3. The planar limit of this free energy can be calculated 
at weak 't Hooft coupling by using perturbation theory. If the superconformal field theory has 
an AdS dual, its value at strong coupling should be given by the regularized gravity action in 
AdS<4. In the case of ABJM theory one obtains in this way 



is the 't Hooft parameter of ABJM theory. In [32] it was shown, from the matrix model represen- 
tation of [61], that the planar free energy of ABJM theory on the three-sphere can be computed 
explicitly at all values of the 't Hooft coupling, and it gives a non-trivial interpolating function 
between the perturbative, weak coupling result, and the strong coupling result quoted in (1.1), 
confirming in this way the prediction of the large N AdS dual. This gives a new, powerful 
test of the AdS/CFT correspondence which has been generalized to many other superconformal 
Chern-Simons-matter theories with AdS duals. 

The goal of these lectures is to provide a simple derivation of the exact planar free energy of 
ABJM theory, first presented in [32] . To begin with, we explain how to obtain the weak coupling 
result in field theory and the strong coupling result in AdS supergravity in (1.1). Then we show 




(1.1) 
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N 
~k 



(1.2) 
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how the recent progress in [61, 32] makes possible to obtain the exact function of the 't Hooft 
coupling interpolating between these two results. To do this, we present the localization technique 
of [61] and the matrix model techniques of [47, 75, 32]. The derivation of the interpolating 
function in [32] used results from the special geometry of Calabi-Yau manifolds, but one can 
obtain it by using only standard ideas in large TV matrix models. This is what we do here. 

The purpose of these lectures is mainly pedagogical. We try to provide detailed accounts 
of the background knowledge underlying the calculations in [61, 32]. For example, we explain 
in some detail how to perform perturbative calculations in Chern-Simons theory and how to 
renormalize the gravity action in AdS spaces. We don't provide however a full overview of the 
large N string/gauge duality of [6], which can be found in other reviews like [67]. Also, we don't 
review the very recent developments based on localization and matrix models in supersymmetric 
Chern-Simons-matter theories, but in the last section we summarize some of the most interesting 
results which have been obtained. Given that these developments are still taking place, it is 
probably still premature to cover them appropriately. 

The organization of these lectures is as follows. Section 2 introduces Euclidean supersymmet- 
ric Chern-Simons-matter theories on the three-sphere and their classical properties. In section 3 
we analyze perturbative Chern-Simons theory in some detail and in some generality. In section 
4 we calculate the free energy of Chern-Simons-matter theories on S 3 at one-loop, and we derive 
the weak coupling result for ABJM theory quoted above in (1.1). Next, in section 5, we look at 
the free energy at strong coupling by using the AdS dual, and we explain the basics of holographic 
renormalization of the gravitational action. In section 6 we review the localization computation 
of [61] (incorporating some simplifications in [48]), which leads to a matrix model formulation 
of the free energy of ABJM theory. In section 7 we review some well-established techniques to 
analyze matrix integrals at large N, which are then used to solve the ABJM matrix model in 
section 8. This makes it possible to derive the interpolating function for the planar free energy 
and, in particular, to verify that its strong coupling limit matches the AdS prediction. Finally, 
in section 9 we give a brief summary of some recent developments. An Appendix collects many 
useful facts about harmonic analysis in S 3 which play a key role in the calculations. 



2. Supersymmetric Chern-Simons-matter theories 

In this section we will introduce the basic building blocks of supersymmetric Chern-Simons- 
matter theories. We will work in Euclidean space, and we will put the theories on the three- 
sphere, since we are eventually interested in computing the free energy of the gauge theory in 
this curved space. In this section we will closely follow the presentation of [48]. 

2.1 Conventions 

Our conventions for Euclidean spinors follow essentially [95]. In Euclidean space, the fermions 
ip and i/j are independent and they transform in the same representation of the Lorentz group. 
Their index structure is 

Y> Q , $ a . (2.1) 
We will take 7^ to be the Pauli matrices, which are hermitian, and 

l»v = ^[7^ 1v\ = iW7 P - (2-2) 
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We introduce the usual symplectic product through the antisymmetric matrix 

Cafi = ( ° ° n V (2.3) 



In [95] we have C = — 1 and the matrix is denoted by e a p. The product is 

eA = e a C a p\?. (2.4) 



Notice that 

It is easy to check that 
and in particular 



eyA = ^ 7 (y) 7 a A Q . (2.5) 
eA = Ae, ey^A = -A7^e, (2.6) 



(y»e) A = -eya (2.7) 
We also have the following Fierz identities 

e (e-0) + e (e0) + (ee) V = (2.8) 

and 

e (#) + 2 (ee) $ + (e 7 ^) 7 M e = 0. (2.9) 
2.2 Vector multiplet and supersymmetric Chern Simons theory 

We first start with theories based on vector multiplets. The three dimensional Euclidean N = 2 
vector superfield V has the following content 

V : A„, a, A, A, D, (2.10) 

where is a gauge field, a is an auxiliary scalar field, A, A are two-component complex Dirac 
spinors, and D is an auxiliary scalar. This is just the dimensional reduction of the M = 1 vector 
multiplet in 4 dimensions, and a is the reduction of the fourth component of A^. All fields are 
valued in the Lie algebra q of the gauge group G. For G = U (N) our convention is that q are 
Hermitian matrices. It follows that the gauge covariant derivative is given by 

3 M + i[^,.] (2.11) 

while the gauge field strength is 

= A u - 8 v Ap + A v \. (2.12) 

The transformations of the fields are generated by two independent complex spinors e, e. 
They are given by, 



5Ap = 


^(e 7 ^A- A 7M e), 




5a = 


^(eA-Ae), 




SX = 


-^ u eF^-De + i^eD^ + 


2i „„ 


sx = 


-^eF^ + De-i^Wpa- 




5D = 


- l -e^D,X- l -D,X^e+ l -[eX. 


,cr] + l -[Xe,a 



(2.13) 



i 

6< 
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and we split naturally 

S = 6 e + 6e. (2.14) 

Here we follow the conventions of [48], but we change the sign of the gauge connection: — > 
—A^. The derivative Z? M is covariant with respect to both the gauge field and the spin connec- 
tion. On all the fields, except D, the commutator [<5 e ,<5e] becomes a sum of translation, gauge 
transformation, Lorentz rotation, dilation and R-rotation: 

[S e , Se]A^ =w u d v A^ + i^/i, - D^A, 
[S e , Se]a =ii> M <9 M cr + i[A, a] + pa, 

[S e , 5,]X =iv^X + l -Q^r u X + i[A, A] + ^ P X + aX, 
[6 e , Se]X =w»d^X + ie^y'A + i[A, A] + -p~X - aX, 
[6 e , 5,}D =iv^d^D + i[A, D] + 2pD + ^(crf-fD^e - ef-fD^t), 

where 

=erfe, 

A ^iAfj, + cree, (2.16) 

a^D^efe-efD^e). 

Here, is the spin connection. As a check, let us calculate the commutator acting on a. We 
have, 

[8 e ,6 s ]<T = 8 e (±e\) -Se(-\xe 

= ¥ (~h^ eF ^ ~ De+ ^ ei v) + ^ a (2 17) 



= le^eD^a + pa, 



where we have used (2.7). 

In order for the supersymmetry algebra to close, the last term in the right hand side of 
[5 € , 5e]D must vanish. This is the case if the Killing spinors satisfy 

-f-fD^e = he, -f-fDnDvE = he (2.18) 

for some scalar function h. A sufficient condition for this is to have 

D„e=^e, D,e=^e (2.19) 

and 

h = ~ (2.20) 
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where r is the radius of the three-sphere. This condition is satisfied by one of the Killing spinors 
on the three-sphere (the one which is constant in the left-invariant frame). Notice that, with this 
choice, p in (2.16) vanishes. 

The (Euclidean) SUSY Chern-Simons (CS) action, in flat space, is given by 



Sscs = - J d 3 x Tr (^A A dA + | A 3 - XX + 2Da^J 

= - J d 3 xTr (t^P (a p 8 u A p + ^A p A u A p ^j - ~XX + 2Da^j . 



(2.21) 



Here Tr denotes the trace in the fundamental representation. The part of the action involving 
the gauge connection A is the standard, bosonic CS action in three dimensions. This action was 
first considered from the point of view of QFT in [28], where the total action for a non-abelian 
gauge field was the sum of the standard Yang-Mills action and the CS action. In [97], the CS 
action was considered by itself and shown to lead to a topological gauge theory. 

We can check that the supersymmetric CS action is invariant under the supersymmetry 
generated by S e (the proof for <5g is similar). The supersymmetric variation of the integrand of 
(2.21) is 

(2SA lt d v A p + 2iSA ll A v A p ) e pup - X5X + 2(5D)a + 2D5a 
= -iX^ed v A p e^ p + X ltl eA p A p e pvp 

- X (-ll^F^ -D + i-fD^) e - %\rD p ea ^ 



- i (D^X) yVe + i[Ae, a]a - ^X^D^ea - XeD. 

The terms involving D cancel on the nose. Let us consider the terms involving the gauge field. 
After using (2.2) we find 

^h^F^e = \X lp ee pup d p A v - X 7p ee pup A p A u (2.23) 

which cancel the first two terms in (2.22). Let us now look at the remaining terms. The covariant 
derivative of A is 

D M A = V+^A + iLVA]. (2.24) 
If we integrate by parts the term involving the derivative of A we find in total 

iX^ed^a + \\rd p ea + 1 (-f\) lp e + [A p , X] 7 p ea 
= iX^ed^a + iXj p D p ea + [A p , X]j p ea, 

where we used that 

(VA) 7M e = -h^e. (2.26) 

The derivative of a cancels against the corresponding term in the covariant derivative of a. 
Putting all together, we find 

iA 7 ^ (ZV) <r ~ iA7 M (D p e) a + [A p , \]-fea + \-fe\Ap, a] + i[Ae, a]a. (2.27) 
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The last three terms cancel due to the cyclic property of the trace. This proves the invariance 
of the supersymmetric CS theory. 

In the path integral, the supersymmetric CS action enters in the form 

exp ( ^S SCS ] (2.28) 



\4tt 

where k plays the role of the inverse coupling constant and it is referred to as the level of the CS 
theory. In a consistent quantum theory, k must be an integer [28] . This is due to the fact that the 
Chern-Simons action for the connection A is not invariant under large gauge transformations, 
but changes by an integer times 8ir 2 . The quantization of k guarantees that (2.28) remains 
invariant. 

Of course, there is another Lagrangian for vector multiplets, namely the Yang-Mills La- 
grangian, 



£ym = Tr 



\ F , V F^ + l -D,oD»o + \ (p + °) 2 + \\^D,\ + L\[j, A] - ^AA 



(2.29) 



In the flat space limit r — > oo, this becomes the standard (Euclidean) super Yang-Mills theory in 
three dimensions. The Lagrangian (2.29) is not only invariant under the SUSY transformations 
(2.13), but it can be written as a superderivative, 

ee£ YM = 6 ? <5 e Tr(-AA - 2L>ct) . (2.30) 

This will be important later on. 

2.3 Supersymmetric matter multiplets 

We will now add supersymmetric matter, i.e. a chiral multiplet in a representation R of the 
gauge group. Its components are 

$ : <f>, 4>, V, i>, F, F. (2.31) 
The supersymmetry transformations are 

5(f) =eip, 

2Ai 

Stp =\^eD^ + ieacp + —^D^ecf) + eF, 

Sip =\^eD^ + ifae + ^^fD^e + Fe, (2 ' 32) 
5F =e{yfD^ - W> - SX<t>) + ^(2A - l^efY, 
5F =e(i^D^ - i^a + i<M) + ^(2A - l^eVV, 

where A is the possible anomalous dimension of 4>. For theories with N > 3 supersymmetry, the 
field has the canonical dimension 

A=\, (2.33) 
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but in general this is not the case. 

The commutators of these transformations are given by 

[S e , 5e}4> = iv^d^cp + iA(f) + Ap<f) - Aa(f), 
[6 e , 5S = iifdnt - i4>A + A P 4> + Aa4>, 

[S e , S E ]iP = Wd^ + le^VV + iA^ + (a + pi; + (1 - A)aV, 

[<5 e , ^ = Wdrf + ~ tyA + (a + ^ pj + (A - l)a^, 

[<5 e , <5 E ]F = ii^fyF + iAF + (A + + (2 - A)aF, 

[8 e , S g ]F = PfidpF - iFA + (A + l)pF + (A - 2)aF. (2.34) 

The lowest components of the superfields are assigned the dimension A and R-charge =fA. The 
supersymmetry algebra closes off-shell when the Killing spinors e, e satisfy (2.18) and h is given 
by (2.20). As a check, we compute 

/ 2iA \ 2iA (2-35) 

= e (i^eD^ + ieacj) + — ^ {D^e) = iu"!^ + iaee + — (e 7 M ^V) - 

which is the wished- for result. 

Let us now consider supersymmetric Lagrangians for the matter hypermultiplet. If the fields 
have their canonical dimensions, the Lagrangian 

3__ 

C = D^D^cj) - i^yD^V + 4^2 00 + + i^A^ - i(/>\ip + + + FF (2.36) 

is invariant under supersymmetry if the Killing spinors e, e satisfy (2.18), with /i given in (2.20). 
The quadratic part of the Lagrangian for <j> gives indeed the standard conformal coupling for a 
scalar field. We recall that the action for a massless scalar field in a curved space of n dimensions 
contains a coupling to the curvature R given by 

S = J d n x^g~ {g^d^d^ + , (2.37) 

where £ is a constant. This action is conformally invariant when 

If the spacetime is an n-sphere of radius r, the curvature is 

*=?*£i>, (2.39) 



and the conformal coupling of the scalar leads to an effective mass term of the form 

n(n-2) 2 

4r 2 ^ 



(2.40) 
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which in n = 3 dimensions gives the quadratic term for (j) in (2.36). 
If the fields have non-canonical dimensions, the Lagrangian 

Anat =DjDl X 4> + 0CT 2 (/> + 1 ( 2A - 1 ) ^ + A ( 2 ~ + ^£,0 + ^ F 

r r (2.41) 

- 2A - 1 - - — 

— iijj'y^ D nip + iipcrip V'V' + i^A</> — i^A^? 

2r 

is supersymmetric, provided the parameters e,e satisfy the Killing spinor conditions (2.19). The 
Lagrangian (2.41) is not only invariant under the supersymmetries 5 6y e, but it can be written as 
a total superderivative, 

ee £ mat = 5,5 e - 2i0<r</> + 2(A ~ . (2.42) 



2.4 ABJM theory 

The theory proposed in [6, 5] to describe N M2 
branes is a particular example of a supersymmetric 
Chern-Simons theory. It consist of two copies of 
Chern-Simons theory with gauge groups U{N\), 
U(N2), and opposite levels k, —k. In addition, we 
have four matter supermultiplets $j, i = 1, • • • , 4, 
in the bifundamental representation of the gauge 
group U(Ni) x U{N2). This theory can be repre- 
sented as a quiver, with two nodes representing the 
Chern-Simons theories, and four edges between 
the nodes representing the matter supermultiplets, 
see Fig. 1. In addition, there is a superpotential 
involving the matter fields, which after integrat- 
ing out the auxiliary fields in the Chern-Simons- 
matter system, reads (on M 3 ) 




$;=i,... ,4 

Figure 1: The quiver for ABJ(M) theory. The 
two nodes represent the U(Ni t 2) Chern-Simons 
theories (with opposite levels) and the arrows 
between the nodes represent the four matter 
multiplets in the bifundamental representation. 



W = — Tr ($!$t$ 3 $T _ $ 1$ T $3$ T j 5 (2 4 3) 

where we have used the standard superspace notation for M = 1 supermultiplets [95]. 
3. A brief review of Chern-Simons theory 

Since one crucial ingredient in the theories we are considering is Chern-Simons theory, we review 
here some results concerning the perturbative structure of this theory on general three-manifolds. 
These results were first obtained in the seminal paper by Witten [97] and then extended and 
refined in various directions in [41, 58, 87, 88, 1,2]. Chern-Simons perturbation theory on general 
three-manifolds is a important subject in itself, hence we will try to give a general presentation 
which might be useful in other contexts. This will require a rather formal development, and 
the reader interested in the result for the one-loop contribution might want to skip some of the 
derivations in the next two subsections. 
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3.1 Perturbative approach 

In this section, we will denote the bosonic Chern-Simons action by 

S = --!-[ Tr(A AcL4+ -A A A A A) (3.1) 
4vr i M \ 3 / 

where we use the conventions appropriate for Hermitian connections, and we included the factor 
1 /4-7T in the action for notational convenience. The group of gauge transformations Q acts on the 
gauge connections as follows, 

A -»• A u = UAU- 1 - iUdU- 1 , UeQ. (3.2) 

We will assume that the theory is defined on a compact three- manifold M. The partition function 
is defined as 

Z{M) = / PA]e ikS (3.3) 

where we recall that k G Z. 

There are many different approaches to the calculation of (3.3), but the obvious strategy 
is to use perturbation theory. Notice that, since the theory is defined on a compact manifold, 
there are no IR divergences and we just have to deal with UV divergences, as in standard QFT. 
Once these are treated appropriately, the partition function (3.3) is a well-defined observable. 
In perturbation theory we evaluate (3.3) by expanding around saddle-points. These are flat 
connections, which are in one-to-one correspondence with group homomorphisms 

7Ti(M) -> G (3.4) 

modulo conjugation. For example, if M = S 3 /Z p is the lens space L(p, 1), one has iri(L(p, 1)) = 
Z p , and flat connections are labelled by homomorphisms Z p — > G. Let us assume that these are 
a discrete set of points (this happens, for example, if M is a rational homology sphere, since in 
that case tt\ (M) is a finite group) . We will label the flat connections with an index c, and a flat 
connection will be denoted by A^ . Each flat connection leads to a covariant derivative 

d Aic) =d + iL4( c V], (3.5) 

and flatness implies that 

d 2 A(c) =iF A(c) =0. (3.6) 
Therefore, the covariant derivative leads to a complex 

-> fi°(M,g) ^ ^(M,g) ^ n\M,g) ^ n 3 (M,g). (3.7) 

The first two terms in this complex have a natural interpretation in the context of gauge theories: 
fi°(M, g) is the Lie algebra of the group of gauge transformations, and we can write a gauge 
transformation as 

U = (M, (f>en°(M,g). (3.8) 
The elements of J7°(M, g) generate infinitesimal gauge transformations, 

5 A = -d A (f>. (3.9) 
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Ker d A ( 



c) 



Imd A(c) 



d 



AM 



n°(M, g ) 



AM 



Imd 



AM 



Kerd 1 ", 



AM 



n\M, s ) 



Figure 2: The standard elliptic decomposition of f2°(M, g) and ®}{M, g). 

The second term, f2 1 (M, g), can be identified with the tangent space to the space of gauge 
connections. The first map in the complex (3.7) is interpreted as (minus) an infinitesimal gauge 
transformation in the background of A^ c \ 

We recall that the space of g- valued forms on M has a natural inner product given by 



(a,b) 



M 



Tr (a A *6), 



(3.10) 



where * is the Hodge operator. With respect to this product we can define an adjoint operator 
on g-valued p-forms in the same way that is done for the usual de Rham operator, 



d^ (c) = (-l) 3(1+p)+1 *d A(c) *. 
We then have the orthogonal decompositions (see Fig. 2) 

n°(M,g) =Kerd A(c) elmd^, 

fi^g) =Kerd^ (c) ®Imd #) . 
These decompositions are easily proved. For the first one, for example, we just note that 



a G Kerd A(c) => (d A ( c )a, <j>) = (a, d 



A( c > 



0, 



therefore 



(Kerd^c))" 1 = Imd^ (c) . 
One also has the analogue of the Laplace-Beltrami operator acting on p-forms 

A #) = d \(c ) d A (c ) +d A(c) d\ (r . ) . 

In the following we will assume that 

H 1 (M,g) = 0. 



(3.11) 

(3.12) 

(3.13) 
(3.14) 

(3.15) 
(3.16) 



This means that the connection is isolated. However, we will consider the possibility that 
has a non-trivial isotropy group H c . We recall that the isotropy group of a connection A^ 
is the subgroup of gauge transformations which leave A^ invariant, 

H c = {0 e g\HA (c) ) = A< c >}. (3.17) 

The Lie algebra of this group is given by zero- forms annihilated by the covariant derivative (3.5), 

Lie(^ c ) = H°(M, g) = Ker d A(c) , (3.18) 
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which is in general non-trivial. A connection is irreducible if its isotropy group is equal to the 
center of the group. In particular, if is irreducible one has 

tf°(M,g) = 0. (3.19) 

It can be shown that the isotropy group % c consists of constant gauge transformations that leave 

A (c) 

invariant, 

4>A^4>- 1 = A^\ (3.20) 

They are in one-to-one correspondence with a subgroup of G which we will denote by H c . 

In the semiclassical approximation, Z(M) is written as a sum of terms associated to saddle- 
points: 

Z{M) = ^Z {c \M), (3.21) 

c 

where c labels the different flat connections A^ on M. Each of the (M) will be a perturbative 
series in l/k of the form 

f oo ^ 



Z^\M) = z[% oop (M)e W \J2sr } k- e , (3.22) 



where is the (I + l)-loop contribution around the flat connection A^ c \ In order to derive 
this expansion, we split the connection into a "background", which is the flat connection A^ c \ 
plus a "fluctuation" B: 

A = A ic) + B. (3.23) 

Expanding around this, we find 

S(A) = S(A®) + S(B), (3.24) 

where 

S(B) = -^f Tr(B/\d A(c) B + ^B 3 ). (3.25) 
4vr J M 3 

The first term in (3.24) is the classical Chern-Simons invariant of the connection A^ c \ Since 
Chern-Simons theory is a gauge theory, in order to proceed we have to fix the gauge. We will 
follow the detailed analysis of [1]. Our gauge choice will be the standard, covariant gauge, 

g A(c) (B) = d\ (c) B = (3.26) 

where g A ( c ) is the gauge fixing function. We recall that in the standard Fadeev-Popov (FP) 
gauge fixing one first defines 

A aI = / VUS (9a(0 (B U )) , (3.27) 
and then inserts into the path integral 



VU5{g A(c) (B u )) 



A A( c) (B) . (3.28) 
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The key new ingredient here is the presence of a non-trivial isotropy group T-L c for the flat 
connection . When there is a non-trivial isotropy group, the gauge-fixing condition does not 
fix the gauge completely, since 



g A{c) (B<t>) = 4>g A(c) (B)(t>- 1 , <f>€H c , 



(3.29) 



i.e. the basic assumption that g(A) = only cuts the gauge orbit once is not true, and there is 
a residual symmetry given by the isotropy group. Another way to see this is that the standard 
FP determinant vanishes due to zero modes. In fact, the standard calculation of (3.27) (which 
is valid if the isotropy group of is trivial) gives 



det 



Sg A (o)(B u ) 



5U 



-i 



det d^ (c) dA 



(3.30) 



But when % c ^ 0, the operator d A ( c ) has zero modes due to the nonvanishing of (3.18), and 
the FP procedure is ill-defined. The correct way to proceed in the calculation of (3.27) is to 
split the integration over the gauge group into two pieces. The first piece is the integration 
over the isotropy group. Due to (3.29), the integrand does not depend on it, and we obtain a 
factor of Vol("H c ). The second piece gives an integration over the remaining part of the gauge 
transformations, which has as its Lie algebra 



(Kerd A(c) )- 



(3.31) 



The integration over this piece leads to the standard FP determinant (3.30) but with the zero 
modes removed. We then find, 



A A l c) (B) = Vol(W c ) 



detd^ (c) dA 



(Kerd (c) ) 



(3.32) 



This phenomenon was first observed by Rozansky in [87], and developed in this language in [1]. 
As usual, the determinant appearing here can be written as a path integral over ghost fields, 
with action 

S ghosts (C,C,B) = (C,d\ (c) d A C), (3.33) 
where C, C are Grassmannian fields taking values in 

(Kerd^J-K (3.34) 

The action for the ghosts can be divided into a kinetic term plus an interaction term between 
the ghost fields and the fluctuation B: 



S ghosts (C,C,B) = (C,A A(c) C)+i(C,d Aic) [B,C}}. 
The modified FP gauge-fixing leads then to the path integral 
z (c) (M) = e ifcS(A«0) f VBe ikS ^A A(c) (B)S (d\ M B) 



(3.35) 



Vol(Hc) A}i(M,g) 



f VB5(d\, c) B) [ T>CVCe ikS W- s z h °^ c > c > B \ 



(3.36) 
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Finally, we analyze the delta constraint on B. Due to the decomposition of J7 1 (M, g) in (3.12), 
we can write 

B = d A(c)( j) + B', (3.37) 

where 

0e(Kerd A(c) )\ B' eKerd^ (c) . (3.38) 

The presence of the operator d A ( c ) in the change of variables (3.37) leads to a non-trivial Jacobian. 
Indeed, we have 

llBfH^A^ + HB'H 2 , (3.39) 
and the measure in the functional integral becomes 

VB = (det' A° (c) ) 3 V(f)VB', (3.40) 

where the ' indicates, as usual, that we are removing zero modes. Notice that the operator in 
(3.40) is positive-definite, so the square root of its determinant is well-defined. We also have 
that 

5(d\ (c) B) =5{A% c) 0) = (det' A^))" 1 5(0), (3.41) 
which is a straightforward generalization of the standard formula 

S(ax) = j^S(x). ( 3 - 42 ) 

We conclude that the delta function, together with the Jacobian in (3.40), lead to the the following 
factor in the path integral: 

(det'A° (c) p. (3.43) 

In addition, the delta function sets (p = 0. It only remains the integration over B' , which we 
relabel B' — >■ B. The final result for the gauge-fixed path integral is then 



p ifcS(A( c )) iff - 

Z< C )(M) = \—— (det' A° (e) )~ 5 / VB VCV c^s{B)-s, hosts {c,c,B)_ 

(3.44) 

This is starting point to perform gauge-fixed perturbation theory in Chern-Simons theory. 
3.2 The one loop contribution 

We now consider the one-loop contribution of a saddle-point to the path integral. This has been 
studied in many papers [41, 58, 87, 88]. We will follow the detailed presentation in [2]. 

Before proceeding, we should specify what is the regularization method that we will use to 
define the functional determinants appearing in our calculation. A natural and useful regular- 
ization for quantum field theories in curved space is zeta- function regularization [50]. We recall 
that the zeta function of a self-adjoint operator T with eigenvalues \ n > is defined as 

Ct(*) = E A « S - ( 3 - 45 ) 

n 

Under appropriate conditions, this defines a meromorphic function on the complex s-plane which 
is regular at s = 0. Since 

-Ct(0) = ^logA n (3.46) 
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we can define the determinant of T as 



det(T) = e^T 



&(0) 



(3.47) 



This is the regularization we will adopt in the following. It has the added advantage that, 
when used in Chern-Simons theory, it leads to natural topological invariants like the Ray-Singer 
torsion, as we will see. 

The main ingredients in the one-loop contribution to the path integral (3.44) are the deter- 
minants of the operators appearing in the kinetic terms for B, C and C. Putting together the 
determinant (3.43) and the determinant coming from the ghost fields, we obtain 



A 1/2 



(det'A° (c) ) 



(3.48) 



since the ghosts are restricted to (3.34) and their determinant is also primed. The operator 
appearing in the kinetic term for B is iQ/2, where 



^ = ~2tt * dj4(c) 

is a self-adjoint operator acting on fi 1 (M, g) which has to be restricted to 

Kerd A( C ) = ( Imd A(=)) ± 
due to the gauge fixing. Notice that, if (3.16) holds, one has 

fl^A^g) = =► Kerd A(c) = Imd^c), 



(3.49) 



(3.50) 



(3.51) 



and due to the restriction to (3.50), Q has no zero modes. However, the operator Q is not 
positive definite, and one has to be careful in order to define its determinant. We will now do 
this following the discussion in [97, 11]. A natural definition takes as its starting point the trivial 
Gaussian integral 



f°° ( 
I dx exp I i 

J — oo \ 




ITT 



exp —sign A 



(3.52) 



If we want to have a natural generalization of this, the integration over B should be 

-1/2 



ITT 



exp ( — sign(Q) 



det 



Q_ 

2tt 



(3.53) 



In order to compute the determinant in absolute value, we can consider the square of the operator 
— * d^( c ) , which is given by 



*d j4(c) * d A ( c 



(3.54) 



and it is positive definite when restricted to (3.50). It is the Laplacian on one- forms, restricted 
to (3.50). We then define 



det 



Q 

'2vr 



det' 



(3.55) 
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where we have subtracted the zero-modes (i.e., we consider the restriction to (3.50)). In order 
to take into account the sign in (3.53), we need the 77 invariant of the operator — * d^( c ). We 
recall that the r\ invariant is defined as 

vt(s) = 7 \s - ( 1 _y ( 3 - 56 ) 

i ( X j ) ^ [ \) 

where are the strictly positive (negative, respectively) eigenvalues of T. The regularized 
difference of eigenvalues is then t]t(0). In our case, this gives 

rK^ (c) )^-*d A(c) (0). (3.57) 



Finally, we have to take into account that for each eigenvalue of the operator (3.54) we have a 
factor of 

-1/2 

(3.58) 



(-) 
Uvr 2 J 



appearing in the final answer. The regularized number of eigenvalues of the operator is simply 
restricted again to (3.50). Putting all together we obtain, 



^i)" = (i) " C(A(C)V2 exp(^(^ (c) )j (det'd^d^)-. (3.60) 

Here we assumed that k > 0. For a negative level —k < the answer is still given by (3.60), but 
the phase involving the eta invariant has the opposite sign. We can now combine this result with 
(3.48). The quotient of the determinants of the Laplacians gives the square root of the so-called 
Ray-Singer torsion of the flat connection [85], 



(det'A° (e) )* 
(det' d^ {c) d A ( c) ) 



4 



'(M,A( C )). (3.61) 



This was first observed by Schwarz in the Abelian theory [89]. When the connection is 
isolated and irreducible, this quotient is a topological invariant of M, but in general it is not. 
However, for a reducible and isolated flat connection, the dependence on the metric is just given 
by an overall factor, equal to the volume of the manifold M: 

t' r {M, A^) = (vol(M)) dim(Wc) t r (M, A^). (3.62) 

where tr(M, A^) is now metric- indepedendent. For an explanation of this fact, see for example 
Appendix B in [42]. However, this volume, which is a metric-dependent factor, cancels in the 
final answer for the one-loop path integral [2]. The isotropy group T-L c is the space of constant 
zero forms, taking values in a subgroup H c C G of the gauge group. Each generator of its Lie 
algebra has a norm given by its norm as an element of g, times 

^ 1 = (vol(M)) 1/2 . (3.63) 
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Therefore, 

vol(W c ) = (vol(M)) dim(Wc)/2 vo\(H c ), (3.64) 

and 

VOl(W c ) VOl(fl" c ) 1 ■ J 

which does not depend on the metric of M. Finally, in order to write down the answer, we take 
into account that, for isolated flat connections, ((A^) can be evaluated as [3] 

C(^l (c) ) = dim#°(M,g). (3.66) 

Putting everything together, we find for the one-loop contribution to the path integral 

Z [\ op (M) = (^) ^'H* ^ (M ,^)). (3.67) 

As noticed above, this expression is valid for k > 0. For a negative level —k < 0, the phase 
involving the gravitational rj invariant has the opposite sign. It was pointed out in [97] that this 
phase can be written in a more suggestive form by using the Atiyah-Patodi-Singer theorem, 
which says that 

11 {A^) = t 1 {Q) + ^S{A^). (3.68) 

Here y is the dual Coxeter number of G (for U(N), y = N), and r?(0) is the eta invariant of the 
trivial connection. Let us denote by 

d G = dim(G), (3.69) 

the dimension of the gauge group. The operator involved in the calculation of r/(0) is just do 
copies of the "gravitational" operator — * d, which is only coupled to the metric. We can then 
write 

77(0) = d GVgrav , (3.70) 
where r/ grav is the "gravitational" eta invariant of — * d . We then find, 



1 / h \ -|dimff°(M,g) . 

Z Soop( M ) = ( 4^2 ) ^W^^yJrnMAM). (3.71) 



-|dim H°(M,g) 

VOl(H c ) [4^ J 

This formula exhibits a one-loop renormalization of k 

k^k + y (3.72) 

which is simply a shift by the dual Coxeter number [97]. However, different regularizations of 
the theory seem to lead to different shifts [10]. 

When = is the trivial flat connection, one has that H c = G, where G is the gauge 
group, and the cohomology twisted by A^ reduces to the ordinary cohomology. The Ray-Singer 
torsion is the torsion tr(M) of the ordinary de Rham differential, to the power da- We then 
obtain, for the trivial connection, 

^i-ioop(M) = ^ J e?**" (r R (M)) d ^ 2 . (3.73) 
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As explained in detail in [97], the phase in (3.71) and (3.73) involving the r? invariant is 
metric-dependent. In constructing a topological field theory out of Chern-Simons gauge theory, 
as in [97], one wants to preserve topological invariance, and an appropriate counterterm has 
to be added to the action. The counterterm involves the gravitational Chern-Simons action 
S(oj), where oj is the spin connection. However, this action is ambiguous, and it depends on a 
trivialization of the tangent bundle to M. Such a choice of trivialization is called a framing of 
the three-manifold. The difference between two trivializations can be encoded in an integer s, 
and when one changes the trivialization the gravitational Chern-Simons action changes as 

S(u) -> S(u) + 2vrs, (3.74) 

similarly to the gauge Chern-Simons action. According to the Atiyah-Patodi-Singer theorem, 
the combination 

r*- + j^r (3 - 75) 

is a topological invariant. It depends on the choice of framing of M, but not on its metric. There- 
fore, if we modify (3.73) by including in the phase an appropriate multiple of the gravitational 
Chern-Simons action, 



exp I ^d G ry grav ] ->■ exp 



Ma [ + 1 S(UJ) 



12 2vr 



(3.76) 



the resulting one-loop partition function is a topological invariant of the framed three-manifold 
M. If we change the framing of M by s units, the above factor induces a change in the partition 
function which at one-loop is of the form 

Z(M) ^exp^nis—^j Z(M). (3.77) 

One of the most beautiful results of [97] is that Chern-Simons theory is exactly solvable on 
any three-manifold M, and its partition function can be computed exactly as a function of k, for 
any gauge group G, by using current algebra in two dimensions. In particular, one can compute 
the exact change of the partition function under a change of framing, and one finds 

Z(M) -> exp (2™ • ^) Z(M), (3.78) 

where 

c=P^. (3.79) 
k + y 

A detailed explanation of the exact solution of CS theory would take us too far, and we refer the 
reader to the original paper [97] or to the presentation in [73] for more details. We will however 
list later on the relevant results when M = S 3 . 



4. The free energy at weak coupling 

The partition function of a CFT on § 3 should encode information about the number of degrees 
of freedom of the theory, in the sense that at weak coupling it should scale as the number N 
of elementary constituents. This follows simply from the factorization property of the partition 
function in the absence of interactions: 

Z(8 3 ,M) ~ {Z{«&\l)) N . (4.1) 
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For example, a gauge theory with gauge group U(N) has at weak coupling TV 2 degrees of freedom, 
and we should expect the free energy on the three-sphere to scale in this regime as 



F(S 3 ) ~ 0{N 2 ). (4.2) 

Of course, at strong coupling this is not necessarily the case. 

In this section we will compute the partition function on § 3 of supersymmetric Chern- 
Simons-matter theories in the weak coupling approximation, i.e. at one loop. First, we will do 
the computation in Chern-Simons theory, and then we will consider the much simpler case of 
supersymmetric matter multiplets. 



4.1 Chern— Simons theory on S 3 

In the previous section we presented the general procedure to calculate the one-loop contribu- 
tion of Chern-Simons theory on any three-manifold, around an isolated flat connection. This 
procedure can be made very concrete when the manifold is S 3 . In this case there is only one 
flat connection, the trivial one = 0, and we can use (3.73). Therefore, we just have to 
compute the Ray-Singer torsion r(S 3 ) for the standard de Rham differential, i.e. the quotient of 
determinants appearing in (3.61) with = (a similar calculation was made in Appendix A 
of [42]). 

We endow S 3 with its standard metric (the one induced by its standard embedding in M 4 
with Euclidean metric), and we choose the radius r = 1 (it is easy to verify explicitly that the 
final result is independent of r). The determinant of the scalar Laplacian on the sphere can be 
computed very explicitly, since its eigenvalues are known to be (see the Appendix) 

X n = n(n + 2), n = 0,l,--- (4.3) 

where n is related to j in (A. 45) by n = 2j. The degeneracy of this eigenvalue is 

d n = (n+l) 2 . (4.4) 

Removing the zero eigenvalue just means that we remove n = from the spectrum. To calculate 
the determinant we must calculate the zeta function, 

°° d °° (n+1) 2 °° m 2 

Cao(s) = m =y £ wn +2) r = ^ 2 w=w (45) 

n=l " 71=1 v v m=2 x ' 

This zeta function can not be written in closed form, but its derivative at s = is easy to 
calculate. The calculation can be done in many ways, and general results for the determinant 
of Laplacians on S m can be found in for example [84, 96]. We will follow a simple procedure 
inspired by [80]. We split 

m 2 1 s 



(m 2 - l) s m 2 ^- 1 ) 



jS + — s + R ( m ' s )> ( 4 - 6 ) 



where 

m 2 I 



R( m > «) = 7 — 2 577-TT 2~s ( 47 ) 

(m z — I) m 2 y s l > m zs 



-19- 



which decreases as m~ 2s ~ 2 for large m, and therefore leads to a convergent series for all s > 
—1/2 which is moreover uniformly convergent. Therefore, it is possible to exchange sums with 
derivatives. The derivative of R(m, s) at s = can be calculated as 



dR(m, s) 



ds 



= —1 — m log ( 1 — 



s=0 



777/ 



(4.8) 



The sum of this series can be explicitly calculated by using the Hurwitz zeta function, and one 
finds 



-E 

m=2 



1 + m log ( 1 



1 



771/ 



- fog(vT). 



We then obtain 



m=2 



Cao(s) = C(2* - 2) - 1 + s (C(2s) - 1) + Y, R ( m > s )> 

where ((s) is Riemann's zeta function, and 

-Clo(0) = log(7T)-2C / (-2). 
The final result can be expressed in terms of C(3), since 

C(3) 



C("2) 



4vr 2 



We conclude that the determinant of the scalar Laplacian on S 3 is given by 

C(3) 



log det' A u = log(vr) + 



2tt 2 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



We now compute the determinant in the denominator of (3.61). We must consider the space 
of one- forms on S 3 , and restrict to the ones that are not in the image of d. These forms are 
precisely the vector spherical harmonics, whose properties are reviewed in the Appendix. The 
eigenvalues of the operator d^d are given in (A. 52), and they read 



A n = (n + 1) 2 , n = l,2,- 



(4.14) 



with degeneracies 



d n = 2n(n + 2). (4.15) 
The zeta function associated to this Laplacian (restricted to the vector spherical harmonics) is 

77l 2 — 1 



= E = 2 E ^ = 2 C(2* - 2) - 2C(2,), (4.16) 

n=l \ ' ' m=l 



and 



log det' A 1 = -4C'(-2) - 21og(27T) = -21og(27r) + 



C(3) 



7T^ 



Here we have used that 



C'(o) 



log(27T) 



(4.17) 
(4.18) 
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We conclude that 



logr^S 3 ) = log det' A^ ) - - log det' = log(2^ 2 ). (4.19) 



This is in agreement with the calculation of the analytic torsion for general spheres in for example 
[94]. In view of (3.62), and since 

vol(§ 3 ) = 2vr 2 , (4.20) 

we find 

tr(S 3 ) = 1. (4.21) 
One can also calculate the invariant (3.59) directly, since this is do times 

C A i(0) = -2C(0) = l, (4.22) 

and it agrees with (3.66). Since the eigenvalues of *d on the vector spherical harmonics come in 
pairs with the same absolute value but opposite signs (see (A.51)), ry grav = 0. We conclude that, 
for k > 0, 

d G 

, 3 > l / * \~ 



In particular, for G = U(N) we have 

1 / k 



z >-'^ = ^uW){^) (4 ' 24) 

The volume of U(N) is given by 

voKu m = (4 . 25) 

where (72(2) is the Barnes function, defined by 

G 2 (z + l) = T(z)G 2 (z), G 2 (l) = l. (4.26) 

Notice that 

G 2 (N + 1) = (N - 1)\(N - 2)! • • • 2!1!. (4.27) 

As we mentioned before, the partition function of Chern-Simons theory on S 3 can be com- 
puted exactly for any gauge group. In the case of U(N), the answer is, for k > [97] 

1 N ~ 1 / ■ \ N-j 

^(S 3 ) = 7 -^m^ TT I 2sin-^— ] . (4.28) 

y ' (k + N) N / 2 y V k + N ) v ' 

Here an explicit choice of framing has been made, but one can compute the partition function 
for any choice of framing by simply using (3.78). The expansion for large k should reproduce 
the perturbative result, and in particular the leading term should agree with the result (4.24). 
Indeed, we have for k large 



Z(S 3 ) « k~ N / 2 J] (^J =(27t)^ N - 1 H- n2 / 2 G 2 (N + 1), (4.29) 
which is exactly (4.24). 
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4.2 Matter fields 



The supersymmetric multiplet contains a conformally coupled complex scalar and a fermion, 
both in a representation R of the gauge group. The partition function at one loop is just given 
by the quotient of functional determinants 



matter fdet(-i#)V« 



J l— loop 



det A r 



(4.30) 



where cLr is the dimension of the representation, and 

A C = A° + ^ 



(4.31) 



is the conformal Laplacian (we have set again r = 1). We now compute these determinants. 
The eigenvalues of the conformal Laplacian are simply 



n(n + 2) + -, n = 0, l,-- - , 
with the same multiplicity as the standard Laplacian, namely (n + l) 2 . We then have 



(n + If 



E 



m 



-(n(n + 2) + f) s ^{m?-\y 
As in the case of the standard Laplacian, we split 

1 



m 



+ 



where 



( m 2 _ i) s m 2(s-i) 4m 2s 
i? c (m, s) 



+ R c (m,s), 



m 1 s 

(m 2 - ~ m 2 ^" 1 ) ~ 4m 2 *' 



The derivative of R c (m, s) at s = is 

dR c (m, s 



ds 



= — - — m 2 log ( 1 — 



s=0 



1 



The sum of this series can be explicitly calculated as 



-E 



m=l 



1 



+ m log ( 1 — 



4m 2 



4m 2 / 



1 1, , n . 7C(3) 



We then find 



Ca c (s) = C(2s - 2) + Jc(2s) + ^ i2c(m, s), 



171=1 



and we conclude that the determinant of the conformal Laplacian on S 3 is given by 

log det A c = -Ck c (0) = \\og(2) - 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



(4.38) 



(4.39) 
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This is in agreement with the result quoted in the Erratum to [30] 1 . 

Let us now consider the determinant (in absolute value) for the spinor field. We have, using 
(A.66), 



oo 



n=l v 2 / 



After a small manipulation we can write it as 



Cu9|(a) =2-2-^ 



(O.m. 4- 1 ^-2 



' (2m + l) s - 2 ^ (2m + 1, 

m>l V 7 m>l V J (4-41) 
1 

2 



s 



=2 (r^-i)as-2)- l -(r-i)c(s), 



where we have used that 

Ep^ip-E^-E^^i-^cw. (4.42) 

m>0 v ; n>l fc>l v ; 

The regularized number of negative eigenvalues of this operator is given by 0mi(0)/2 and it 
vanishes, so the determinant of the Dirac operator equals its absolute value. We deduce 

log det {-Up) = -C' lp ,(0) = -g^C(3) - \ log 2. (4.43) 

Combining the conformal scalar determinant with the spinor determinant we obtain, 

log det {-ilp) - log det A c = - ^ log 2. (4.44) 

This can be seen directly at the level of eigenvalues. The quotient of determinants is 

n - n M (4.45, 

mil {m 2 -\) m mi( m -2) 

and its regularization leads directly to the result above (see Appendix A of [32]). We conclude 
that 

Z^oop = 2~ dR l 2 . (4.46) 
4.3 ABJM theory at weak coupling 

We can now calculate the free energy on S 3 of ABJM theory. We will restrict ourselves to 
the "ABJM slice" where the two gauge groups have the same rank, i.e. the theory originally 
considered in [6]. We have two copies of CS theory with gauge group U{N) and opposite levels 
k, —k, together with four chiral multiplets in the bifundamental representation of U{N) x U{N). 
Keeping the first term (one-loop) in perturbation theory we find, at one-loop, 

Fabjm^ 3 ) ~ -N 2 log (Jj^j - 2 log {vo\{U{N))) - 2N 2 log(2) (4.47) 



1 Beware: the arXiv version of this paper gives a wrong result for this determinant. 
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where the first two terms come from the CS theories, and the last term comes from the super- 
symmetric matter. Here we assume k > 0. Notice that the theory with opposite level —k gives 
the same contribution as the theory with level k. In order to obtain the planar limit of this 
quantity, we have to expand the volume of U (N) at large N. Using the asymptotic expansion of 
the Barnes function 

N 2 1 3 1 

log G 2 (N + 1) = — log N - — log N - -N 2 + —N log 2ir + C'(-l) 

oo R 

y _ggg 2-2g (44g) 

^ 2 2g{2g-2) 1 > 

where B^g are the Bernoulli numbers, we finally obtain the weakly coupled, planar result 

F A B.m(S 3 ) ~ A^ 2 |log(2vrA) - | - 21og(2)| . (4.49) 



5. Strong coupling and AdS duals 

Some SCFTs in three dimensions have AdS duals given by M-theory /string theory on back- 
grounds of the form 

X = AdS 4 x X 6J , (5.1) 

where Xqj is a six-dimensional or seven-dimensional compactification manifold, depending on 
whether we consider a superstring or an M-theory dual, respectively. One of the consequences of 
the AdS / CFT duality is that the partition function of the Euclidean gauge theory on § 3 should 
be equal to the partition function of the Euclidean version of M-theory/string theory on the dual 
AdS background [98], i.e. 

Zcft(§ 3 ) = Z(X). (5.2) 

In the large N limit, we can compute the r.h.s. in classical (i.e. genus zero) string theory, and 
at strong coupling it is sufficient to consider the supergravity approximation. This means that 
the partition function of the strongly coupled gauge theory on S 3 in the planar limit should be 
given by 

Z CFT (S 3 ) ~ e~ /(AdS4) (5.3) 

where / is the classical gravity action evaluated on the AdS 4 metric. This gives a prediction for 
the strongly coupled behavior of the gauge theory. However, the gravitational action on AdS 4 is 
typically divergent, and it has to be regularized in order to obtain finite results. I will now review 
the method of holographic renormalization and work out two examples: the (related) example of 
the Casimir energy of M = 4 SYM on M x S 3 , and the case of main interest for us, namely the 
free energy of ABJM theory on S 3 . 

5.1 Holographic renormalization 

The gravitational action in an Euclidean space with boundary has two contributions. The first 
contribution is the bulk term, given by the Einstein-Hilbert action 

/bulk = - T7T~n~ I dH+lx ^G(R~ 2A ) (5-4) 
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where Gat is Newton's constant in n + 1 dimensions and G is the (n + l)-dimensional metric. 
The second contribution is the surface term [44] 



^surf 



8ttG n 



dM 



K\-y\V 2 d n x, 



(5.5) 



where dM is the boundary of spacetime, 7 is the metric induced by G on the boundary, and K 
is the extrinsic curvature of the boundary. K satisfies the useful relation (see for example [83]) 



(5.6) 



where n is the normal unit vector to dM, and C n is the Lie derivative along this vector. Both 
actions, when computed on an AdS background, diverge due to the non-compactness of the space. 
For example, after using Einstein's equations, the bulk action of an AdS space of radius L can 
be written as 



n 



'bulk 



/ 



d n+1 xVG 



(5.7) 



8irG N L 2 

which is proportional to the volume of space-time, and it is divergent. 

In order to use the AdS/CFT correspondence, we have to regularize the gravitational action 
in an appropriate way. The procedure which has emerged in studies of the AdS / CFT correspon- 
dence is to introduce a set of universal counterterms, depending only on the induced metric on 
the boundary, which lead to finite values of the gravitational action, energy-momentum tensor, 
etc. This procedure gives values for the gravitational quantities in agreement with the corre- 
sponding quantities computed in the CFT side, and we will adopt it here. It is sometimes called 
"holographic renormalization" and it has been developed in for example [51, 12, 37, 27]. We now 
present the basics of holographic renormalization in AdS. Useful reviews, focused on AdSs, can 
be found in for example [91, 66]. 

An asymptotically AdS metric in n + 1 dimensions with radius L and cosmological constant 



A 



can be written near its boundary at u = as 



n(n — 1) 
2lJ~ 



(5.8) 



ds 2 



+ \gij(u 2 , x)dx l dx j 
u z u z 



The metric gij{u 2 ,x) can be expanded in a power series in u near u = 0, 



9l] (u 2 , x) = sgV) + u 2 gV\x) + u 4 g^(x) + \og{u 2 )h^( 



,2„(2) 



+ 



(5.9) 



(5.10) 



where the first term, gf), is the metric of the CFT on the boundary. The coefficients g\- ^ 

appearing here can be solved recursively in terms of gf^ by plugging (5.10) in Einstein's equations 
One finds, for example [27] 2 



(2) 



n-2 



Rij 



(5.11) 



2 The sign in the curvature is opposite to the conventions in [27], which give a positive curvature to AdS. 
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where Rij and R are the Ricci tensor and curvature of gf^ . The resulting metric is then used to 
compute the gravitational action with a cut-off at u = e which regulates the divergences, 



I e = - 



1 



IQ-kG n 



d n+l x VG [ R + 



M, 



n(n — 1) 
T 2 



f 



dM € 



(5.12) 



Here, M e is the manifold with u > e and a boundary dM t at u = e. To calculate the boundary 
term, we consider the normal vector to the hypersurfaces of constant u, 



u 

T 



(5.13) 



The minus sign is due to the fact that the boundary is at u = 0, so that the normal vector points 
towards the origin. The induced metric is 



L 2 

jijdx t dx :> = —^gij(u 2 ,x)dx l dx-', 



with element of volume 



ii 



s/g. 



(5.14) 
(5.15) 



The intrinsic curvature of the hypersurface at constant u is then 



L 



u 



V9 



nL 



n-1 



U" 



(l - Uid^j y/g. (5.16) 



We then find, 



L = z „ / d n x 



du ^ nL n 1 

Vg- 



8ttG n J " ~ J e u n+1W!J 8ttG N € 
The singularity structure of this regulated action is [51, 27] 

L n-1 



1 ud u yfg 

n 



(5.17) 



1671(7 at 



J d n xJgM (£""0(0) + e- n+2 a {2) + • • • - 2 log(e)a (n) ) + 0(e°). (5.18) 



The logarithmic divergence appears only when n is even. In order to regularize power-type 
divergences in n = 3 and n = 4, it suffices to calculate the first two coefficients, d( ) an d «(2)- 
Let us now calculate these coefficients (the next two are computed in [27]). We first expand, 



so that 



det g = det (l + u 2 Tr (^"V 2 >) +•••), 

VW^x) = yfifi) (l + ^Tr (g^gV) +•••), 
(l - l -ud,)j y^x) = yfi® (l + ^u 2 Tr (g^g®) + • 



(5.19) 



(5.20) 



The regulated Einstein-Hilbert action gives 



nL 



n—l 



87tGtv 



x \ I gy 



1 



+ 



ne n 2(n - 2)e n " 2 



Tr (V°)-V 2 )) 



+ 



(5.21) 



- 26 - 



d n x Jg(°) 



while the regulated Gibbons-Hawking term gives 
nL n - 1 

In total, we find 

L n-1 



1 + ^Tr 



+ 



e WttG n 
and we deduce, 



d n xJg(°) 



2(1 -n) n 2 -5n + A 



(n - 2)e 



n-2 



Tr (^"V 2 )) 



a (o) = 2(1 - n), 

(n-4)(n- 1) 



0(2) = 



n-2 



Tr V 2) ) , n/2. 



(5.22) 



(5.23) 



(5.24) 



The counterterm action is obtained by using a gravitational analogue of the minimal subtraction 
scheme, and it is given by minus the divergent part of I e , 

L n - 1 f .„ /-^T2(n-1) (n-4)(n-l) r 



Ict = 



16nG N 



d n xJg(°) 



(n - 2)e 



+ ■■■ 



+ 21og(e)a (n) 



(5.25) 



As pointed out in [12], we should re- write this in terms of the induced metric in the boundary 
(5.14), evaluated at u = e. From (5.20) we deduce 

,2 



On the other hand, from (5.11) we obtain 



(5.26) 



V y J n-2 V 13 2(n - 1) y 



(0)« 1 = _. 



2(n - 1) 



R 



2(n^lF i?H + 



(5.27) 



where in the first line the Ricci tensor and curvature are computed for gf^ , while in the second line 
the curvature is computed for the induced metric 7. Plugging these results into the counterterm 
action we find 



'ct 



1 + 



4(n- 1) 



Rh} + 



2(n " 1} " 2(n-2) L2fi[7] + 21og(e) °W [7] + • • • 



(5.28) 



- I rPx ^ ( 2 log(£)o '"» l71 + "-ir + 25^2) flH + • • • ) 



2(n-2)" 

which is the result written down in [27, 37] (for Euclidean signature). This is the counterterm 
action which is relevant for AdS in four and five dimensions, and the dots denote higher or- 
der counterterms (in the Riemann tensor of the induced metric) which are needed for higher 
dimensional spaces [12, 37, 27]. The total, regularized gravitational action is then 

I = 4ulk + 4urf + let (5.29) 
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and it yields a finite result by construction. The removal of these IR divergences in the gravita- 
tional theory is dual to the removal of UV divergences in the CFT theory. It can be verified in 
examples that the gravity answers obtained by holographic renormalization match the answers 
obtained in CFT on a curved background after using zeta- function regularization [12]. In the 
next subsection we work out a beautiful example of this matching closely related to the tech- 
niques developed here, namely the Casimir energy for M = 4 super Yang-Mills on M x § 3 , which 
was first derived in [12]. 

5.2 Example 1: Casimir energy in AdSs 

Let us consider an n-dimensional CFT on the manifold S n_1 x SL where the one-circle has 
length /3, and with periodic boundary conditions for the fermions. The supersymmetric partition 
function on this manifold is given by 

^cft (S"- 1 x Sj) = Tr [(-l^e-W 1-1 )] « e"^° (5.30) 

where in the last step we have considered the large (3 limit, and Eq is the energy of the ground 
state, i.e. the Casimir energy on S n_1 . The AdS/CFT correspondence implies that this partition 
function can be obtained by computing the partition function of a superstring/M-theory on a 
manifold of the form (5.1), where now the AdS space has the boundary S™ -1 x §^ [98]. In the 
SUGRA approximation, this can be computed by evaluating the regularized gravity action I 
(5.29), and this should give the planar, strongly coupled limit of the CFT partition function, 

Zcft (S^ 1 x Sj) « e - J ( AdS »+ 1 ). (5.31) 

In order to calculate I, we need an Euclidean AdS metric which is asymptotic to S n_1 x §^ at 
the boundary. The relevant metric is given by [37] 

1 + T2 ) dr 2 + y^j^ + r'dnU (5.32) 

where the boundary is now at r — > oo. The sphere S™" 1 at the boundary has radius L. Here, 
dfi 2 is the metric on the unit n-sphere, and r has length /3. Let us compute the gravitational 
action for this theory, with a cutoff at the boundary dM located at r. We will only include the 
counterterms presented in (5.28), which are sufficient in three and four dimensions. At the end 
of the calculation we will take r — > oo. Let us denote 

V(r) = l + ^. (5.33) 

The normal unit vector to the boundary is 

n = F 1/2 (r)|-, (5.34) 

while the induced metric is 

~/ ij dx i dx j = V{r)dr 2 + r 2 dI2^_!. (5.35) 
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The calculation of the different contributions to the total gravitational action is straightforward: 
8vrGW/ bu i k = -Ir'VoKS"- 1 )/?, 

SvrGWW = -voKS"- 1 )/?^- 1 |^y(r) + \v'{r)\ , (5.36) 
8nG N I ct = voKS"" 1 )^"- 1 + { ^^-) V 1/2 (r). 



The total action is 
8ttG n I 

_ voKS™- 1 )/? 



,71-1 



(n-l)L 2 / r 2 \ . 1N / L 2 \ / L 2 ^ 1/2 ' 



,2 



(5.37) 



Expanding for r — > oo we find, for n = 3, a vanishing action, while for n = 4 (i.e. for AdSs) we 
find 

8 ^«^.|^. ( , 38) 

In this approximation the Casimir energy is then 

Eo = (5.39) 
32Gjv 

Using the standard AdSs/CFT4 dictionary for ftf = 4 super Yang-Mills [71] 

iV2 = (5 - 40) 
we obtain the planar, strong coupling value of the Casimir energy of M = 4 super Yang-Mills on 

Using the technology developed in the previous section it is an easy exercise to compute the 
Casimir energy directly in QFT, at weak coupling. A massless scalar field in a four-dimensional 
curved space with Minkowski signature satisfies the wave equation 

(gVDpDv + &)</> = 0. (5.42) 

For a conformally coupled scalar in 4d, we have 

£ = \- (5-43) 

Let r), x be coordinates for M x S 3 . The metric can be written as 

ds 2 = L 2 (drj 2 - dfll) (5.44) 

where L is the radius of § 3 , and dQ 2 is the element of volume on an § 3 of unit radius. The 
coordinate rj (which is dimensionless) is called the conformal time parameter (see [15], p. 120), 
and it is related to the time coordinate by 

t = Lr). (5.45) 
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We write the wavefunctions in factorized form 

u m (?7,x) = Xm(v)Sm=l{x), (5-46) 
where Sj(x) is a scalar spherical harmonic (A. 41), and it is an eigenfunction of the Laplacian, 

A°S m -i = (m 2 - 1) Sm-i , m = l,2,---. (5.47) 

These functions have degeneracy m 2 . u m satisfies the wave equation (5.42), which after separation 
of variables leads to the following equation for Xm(v) : 

d 2 Xm + (m 2 - 1 + L 2 ^R) Xm = 0. (5.48) 

Here R is the curvature of an S 3 of radius L. In the case of a 4d conformally coupled scalar, the 
above equation reads 

d 2 Xm + m 2 X m = (5.49) 

with the solution 

Xm (v) oc e~ imt / L . (5.50) 
The Casimir energy for this conformally coupled scalar is given by 

^ oo 

e = - m2 • ( m / L ) ( 5 - 51 ) 

m=l 

where the first factor m 2 comes from the degeneracies. The above sum is of course divergent, 
but we can use zeta-function regularization to obtain 

1 00 1 
E ^ = 2L £ m2m ' S = 2L C(S " 2) (5 - 52) 

m=l 

which can be analytically continued to s = — 1. We find in this way 

B _ B( _,)_^ f (_S)_-^. (5.53) 

For Weyl spinors, the Casimir energy is obtained by summing over the modes of the spinor 
spherical harmonics, and with a negative sign due to Fermi statistics, i.e. 

^ oo 

Spinor = Yl 2 < n + + V2), (5-54) 

n=l 

where 2n(n + 1) is the degeneracy of the eigenvalue n + 1/2 of the Dirac operator. This can be 
again regularized by considering 

^spinor(s) = -7)] / n s (5.55) 

L ^1 ( n + 2) 



and analytically continuing it for s = — 1. Using the previous result (4.41) we deduce 
£ spinor (_l) = -I (( 2 - 3 - 1) C (-3) - \ (2- 1 - 1) C(-l)) 17 



960L 



(5.56) 
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Finally, for a gauge field, we have 



oo 



^gau g c = ^^2n(n + 2)(n + l), (5.57) 

71=1 

where n+ 1 is the square root of the energies of the modes (i.e. the square root of the eigenvalues 
(4.14) of the Laplacian, after gauge fixing). This is regularized as 

«.W -\ t -\ t 4 («. - *) - CW) . (-8, 



L """ ■ I. ^ (ll + I ) " / 

m=l 

and one finds 



^ gaug c(-l) = \ (C("3) - C(-l)) = ^ (5-59) 

It follows that the Casimir energy on S 3 x R for a QFT with no conformally coupled real scalars, 
n i/2 Weyl spinors and rt\ vector fields is 

E = ^ (4n + 17n 1/2 + 88m) ■ (5.60) 

This was first shown in [40] by a cutoff regularization of the sum over modes. In the case of 
N = 4 SYM with gauge group U (N) we have 

n = 6iV 2 , m/ 2 = 4iV 2 , m = iV 2 , (5.61) 

and the Casimir energy obtained at weak coupling agrees with (5.41), which is a priori the 
strong-coupling value. This was first observed in [12], and it can be also regarded [27, 90] as a 
consequence of the agreement between the trace anomaly in both sides of the correspondence 
[51]. Indeed, it is known that, for conformal fields in conformally flat backgrounds like S n_1 x S 1 , 
the vev of the energy-momentum tensor is completely determined by the trace anomaly [20], 
which fixes then the value of the Casimir energy on S n_1 [22]. The weak and the strong coupling 
results are the same because these quantities are protected by a non-renormalization theorem 
and do not depend on the 't Hooft coupling A = g 2 N . This will not be the case for the free 
energy on the sphere for 3d SCFTs, which we now compute at strong coupling by using the large 
N AdS dual. 

5.3 Example 2: free energy in AdS4 

We are interested in studying CFTs on S n . Therefore, in the AdS dual we need the Euclidean 
version of the AdS metric with that boundary, which can be written as [37] 

d, 2 = + r 2 d^ 2 (5.62) 

with the notations of the previous subsection. The boundary is again at r — > oo. This metric 
can be also written as [98, 37] 

ds 2 = L 2 (dp 2 + sinh 2 (p)d0 2 ) . (5.63) 

Let us compute the regularized gravitational action (5.29), again with a cutoff at the boundary 
dM located at constant r. The element of volume of the metric G given in (5.62) is 

^ = vmw^ (5 ' 64) 
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where g§n is the metric on an n-sphere of unit radius. The bulk action is just (5.7), i.e. 



-fbulk 



nvol(S r 



dp- 



8ttG n L J '\/ L 2 + p 2- 

To calculate the surface action, we notice that the unit normal vector to dM is 

d 



n = y/l + r 2 /L 2 



dr 



while the induced metric is 



which has the element of volume 



and scalar curvature 



7 = r g E 



We then obtain 



and the surface term is 



^surf 



nr 



n-l 



8itG n 

Finally, the first two counterterms are given by 

Lr n ~ 2 n(n - 1) 



Vl + r 2 /L 2 vol(S n ). 



vol(S n ) 



n-l 



-r" + 



L ' 2(n-2) 
Putting everything together we obtain, 
_ vol(S 



vol(S n ) r n {n - 1) 
8itGn L 



1 + 



n L 2 

2(n-2) r 2 ' 



8ttGnL 



r/L 

nL n I du 

o 



+ r n {n-l) 1 + 



-nr n "Vr 2 + L 2 
n L 2 ~ ~ 



2(n - 2) r 2 

We should now take the limit of this expression when r — >■ oo. For n = 3 we find 



therefore we obtain a /imte action given by 



J = 



vrL 2 



2G 



N 



(5.65) 



(5.66) 

(5.67) 
(5.68) 

(5.69) 

(5.70) 

(5.71) 

(5.72) 



(5.73) 



(5.74) 



(5.75) 



This will give us the strong coupling prediction for the free energy on S 3 of supersymmetric 
Chern-Simons-matter theories with an AdS dual. 
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5.4 ABJM theory and its AdS dual 

In order to compute the free energy of ABJM theory at strong coupling we have to be more precise 
about the gauge/gravity dictionary. We will now write down this dictionary for supersymmetric 
Chern-Simons-matter theories, which was first established for ABJM theory in [6]. We will not 
attempt here to review the derivation of the duality. A pedagogical introduction can be found 
in [67]. 

The AdS duals to the theories we will consider are given by M-theory on 

AdS 4 x X 7 , (5.76) 

where X-j is a seven-dimensional manifold. In the case of ABJM theory, 

X 7 = 8 7 /Z k . (5.77) 

The eleven-dimensional metric and four-form flux are given by the Freund-Rubin background 
(see [36] for a review) 



( ' ls lJ - L X 7 [\ ds - + rls 'A- T I ' 
3.3 



(5.78) 



where WAdS 4 is the volume form with unit radius. The radius Lx 7 is determined by the flux 
quantization condition 

(2tt^) 6 Q= I MiF = 6L% 7 vol(X 7 ). (5.79) 
Jx 7 

In this equation, l p is the eleven-dimensional Planck length. The charge Q is given, at large 
radius, by the number of M2 branes N, but it receives corrections [14, 7]. In ABJM theory we 
have 

Q = N -lAk-jY (5.80) 



24 V fc. 
This extra term comes from the coupling 

C 3 A I s (5.81) 

in M-theory, which contributes to the charge of M2 branes. Here, Ig is proportional to the Euler 
density in eight dimensions, and it satisfies 

(5 - 82) 

where M$ is a compact eight-manifold. In ABJM theory, the relevant eight-manifold is C 4 /Z/%, 
with regularized Euler characteristic 

x(C 4 /Z fe ) =k-\. (5.83) 

This leads to the shift in (5.80). 
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One final ingredient that we will need is Newton's constant in four dimensions. It can 
be obtained from the Einstein-Hilbert action in eleven dimensions, which leads to its four- 
dimensional counterpart by standard Kaluza-Klein reduction, 

-L- /d"xv5TTii„ -> \ ■ J^voKX.) /dSV5«. = f**J&U, (5.84) 

where Gn, Gat denote the eleven-dimensional and the four-dimensional Newton's constant, re- 
spectively, and the volume of Xj is calculated for unit radius. In the resulting Einstein-Hilbert 
action in four dimensions, the metric and scalar curvature refer to an AdS4 space of radius Lx 7 , 
and not Lx 7 /2 as in (5.78). This is the source for the extra factor of 1/4 in the second term. 
Recalling that 

16ttGh = (2ir£ p ) 9 , (5.85) 



we obtain 

1 2V6n 2 Q 3 / 2 1 



G N 9^vol(X 7 ) L 2 Xr ' 
It follows that the regularized gravitational action (5.75) is given by 



(5.86) 



2G N ^ V 27vol(X 
In particular, for ABJM theory we have 



^=03/2/ ^ (5 . 87) 



^fc 1 /2Q3/2_ (5 _ 88) 
o 



In the supergravity and planar approximation we can just set Q = N, and we find indeed that 
the planar free energy is given by 



1 / 2-7T 6 1 

In the case of ABJM theory we finally obtain the strong coupling result stated in (1.1). 



6. Localization 

Localization is an ubiquitous technique in supersymmetric QFT which makes possible to reduce 
an infinite-dimensional path integral to a finite dimensional integral. It features prominently in 
Witten's topological quantum field theories of the cohomological type, where one can argue that 
the semiclassical approximation is exact, see [16, 25] for reviews and a list of references. 

The basic idea of localization is the following. Let 5 be a Grassmann-odd symmetry of a 
theory with action S((f)), where <p denotes the set of fields in the theory. We assume that the 
measure of the path integral is invariant under 5 as well (i.e. S is not anomalous), and that 



S 2 = Cb 



(6.1) 



where Cb is a Grassmann-even symmetry. In a Lorentz- invariant, gauge invariant theory, Cb 
could be a combination of a Lorentz and a gauge transformation. Consider now the perturbed 
partition function 

Z{t) = J V0e- s - t5v , (6.2) 

where V is a Grassmann-odd operator which is invariant under Cb- It is easy to see that Z(t) is 
independent of t, since 

^ = - Jv<t>6V e- s ~ t5v = -Jv<j>8 (Ve- S ~ t5v ) = 0. (6.3) 

Here we have used the fact that 5 2 V = CbV = 0. In the final step we have used the fact that 5 
is a symmetry of the path integral, in order to interpret the integrand as a total derivative. In 
some cases, the integral of the total derivative does not vanish due to boundary terms (a closely 
related example appears in section 11.3 of [78]), but if the integral decays sufficiently fast in field 
space one expects the perturbed partition function Z(t) to be independent of t. This means that 
it can be computed at t = (where one recovers the original partition function) but also for other 
values of t, like t — > oo. In this regime, simplifications typically occur. For example, if 5V has a 
positive definite bosonic part (5V)b, the limit t — > oo localizes the integral to a submanifold of 
field space where 

(5V)b = 0. (6.4) 

It turns out that, in many interesting examples, this submanifold is finite-dimensional. This 
leads to a "collpase" of the path integral to a finite-dimensional integral. It is easy to see that 
this method also makes it possible to calculate the correlation functions of 5-invariant operators. 

In order to see how the method of localization works, let us briefly review a beautiful and 
simple example, namely the field theoretical version of the Poincare-Hopf theorem. 

6.1 A simple example of localization 

The Poincare-Hopf theorem has been worked out from the point of view of supersymmetric 
localization in many references, like for example [16, 25, 70]. Let X be a Riemannian manifold 
of even dimension n, with metric g^ u , vierbein e", and let be a vector field on X. We will 
consider the following "supercoordinates" on the tangent bundle TX 

{^,B„), (6.5) 

where the first doublet represents supercoordinates on the base X, and the second doublet 
represents supercoordinates on the fiber, ip^ and ^ are Grasmann variables. The above super- 
coordinates are related by the Grasmannian symmetry 

6.6 

5^ = 0, 5B„ = 0, 
which squares to zero, S 2 = 0. With these fields we construct the "action" 

S(t) =6*, * = ^ VV (B^ 1 + 2itV + r^VvVV T ) , (6-7) 
and we define the partition function of the theory as 

Zx(t) = j^— [ dzd^dySdSe-SW. (6.8) 
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Using that 

= -T^ - T^g-v, (6.9) 
one finds that, in the resulting theory, is a Gaussian field with mean value 

= -itV - g^T^M" . (6.10) 
If we integrate it out, we obtain an overall factor 

and the action becomes 

^ W " \R p \u^ P ^W - itV^V^r- (6-12) 
We can define orthonormal coordinates on the fiber by using the inverse vierbein, 

Xa = ££YV, (6.13) 

so that the partition function reads 

Z x {t) = 1 J dxdVdxe-T9*.^ ,i ^+3« ofc ^XaX6^^+itv M Ve 5 x-^. ( 6 .i4) 

It is clear that this partition function should be independent of t, since the action can be 
written as 

S(t) = S{0) + tSV, V = l^V. (6.15) 

We can then evaluate it in different regimes: t — > or t — > oo. The calculation when t = is 
very easy, since we just have 

Z X (0) = — L__ ^dxdV,dxei« W^ = ^ J x dxPf(R), (6.16) 

where we have integrated over the Grassmann variables \ a to obtain the Pfaffian of the matrix 
R ab . The resulting top form in the integrand, 

< x ) = (dr^ Pf(jR) (fU7) 

is nothing but the Chern-Weil representative of the Euler class, therefore the evaluation at t = 
produces the Euler characteristic of X, 

Z X (0) = X (X). (6.18) 

Let us now calculate the partition function in the limit t — > oo. We will now assume that 
has isolated, simple zeroes where V^ip^) = 0. These are the saddle-points of the "path 
integral," so we can write Z x (t) as a sum over saddle-points pk, and for each saddle-point we 
have to perform a perturbative expansion. Let be coordinates around the point pk- We have 
the expansion, 

= E ^ • ■■ d ^(Pk)e i ■■■e n - (6.i9) 

n>l 
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After rescaling the variables as 

^r 1/2 V, x^t- 1 /^, (6.20) 

the theory becomes Gaussian in the limit t — > oo, since higher order terms in the fluctuating fields 
€,ip,X contain at least a power t~ 1 / 2 . Interactions are suppressed, and the partition function is 
one- loop exact: 

lim Z x (t) = ^ j^j- 2 I ^d X e-^ k) ^^ + ^^ (6.21) 



t->oo 

Pk 

where we denoted, 

H^ = d a V»\ (6.22) 

'Pk 

Each term in this sum can now be computed as a product of a bosonic Gaussian integral, times 
a Grassmann integral, and we obtain 

1 Het ff(ty 

The equality between (6.18) and (6.23) is the famous Poincare-Hopf theorem. 

Conceptually, the localization analysis in [61] that we will review now is not very different 
from this example, although technically it is more complicated. The key common ingredient 
in the analysis of the t — > oo limit is that the localization locus becomes very simple, and all 
Feynman diagrams involving at least two loops are suppressed by a factor i^ 1 / 2 , so that the 
one- loop approximation is exact. 

6.2 Localization in Chern Simons matter theories: gauge sector 

We are now ready to use the ideas of localization in supersymmetric Chern-Simons-matter 
theories on S 3 , following [61]. The Grassmann-odd symmetry is simply Q, defined by 5 e = eQ, 
where e is the conformal Killing spinor satisfying (2.19). This symmetry satisfies Q 2 = 0, and 
then it is a suitable symmetry for localization. To localize in the gauge sector, we add to the 
CS-matter theory the term 

-tS YM , (6.24) 

which thanks to (2.30) and (2.15) is of the form QV, and its bosonic part is positive definite. 
By the localization argument, the partition function of the theory (as well as the correlators of 
Q-invariant operators) does not depend on t, and we can take t — > oo. This forces the fields to 
take the values that make the bosonic part of (2.29) to vanish. Since this is a sum of positive 
definite terms, they have to vanish separately. We then have the localizing locus, 

F^ = 0, D^a = 0, D+°=0. (6.25) 

The first equation says that the gauge connection must be flat, but since we are on S 3 the 
only flat connection is A^ = 0. Plugging this into the second equation, we obtain 

d^<r = =y- <j = fj , (6.26) 

a constant. Finally, the third equation says that 

D = - — . (6.27) 
r 
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The localizing locus is indeed finite-dimensional: it is just the submanifold where a and D are 
constant Hermitian matrices, and = 0. 

Let us now calculate the path integral over the vector multiplet in the limit t — > oo. We 
have to perform a gauge fixing, and we will choose the standard covariant gauge (3.26) as in the 
case of Chern-Simons theory. The path integral to be calculated is 

— ^(det'A )-! / VA [ DCPCeS s scs-«SYM(4)-s gtote (C,c,A) ) ( 628 ) 

Vbl(G) JKcrdt ./(Kcrd) 1 - 

where C, C are ghosts fields. As in the example of the Poincare-Hopf theorem, we expand the 
fields around the localizing locus, and we set 

1 , 

<t = <7o H — -^a , 
Vt 

D = - a ° + ±_D>, (6.29) 

where the factors of t are chosen to remove the overall factor of t in the Yang-Mills action. In 
the Yang-Mills Lagrangian, only the terms which are quadratic in the fluctuations survive in this 
limit, namely, 

l - I yfg d 3 x Tr T — A^AA^ — [A^, cjq] 2 + d^a'd^a' + (D' + a') 2 

V (6.30) 
+ iA 7 ^V M A + iA[a , A] - ^AA + d^cj , 

where we set r = 1. We are then left with a Gaussian theory, but with non-trivial quadratic 
operators for the fluctuations. In the same way, when we expand (2.21) around the fixed-point 
limit (6.29), we obtain 

t SsCS = ^ Tr K 2 )vol(S 3 ) + <D(t~V 2 ), (6.31) 

so only the first term survives in the t — > oo limit. 

Let us now calculate the path integral when t — > oo. Like in the example of the Poincare- 
Hopf theorem, we just have to compute the one-loop determinants. In this calculation we will 
only take into account the factors which depend explicitly on a®. The remaining, numerical 
factors (which might depend on N, but not on the coupling constant k) can be incorporated 
afterwards by comparing to the weak coupling results. The integral over the fluctuation D' can 
be done immediately. It just eliminates the term (D' + a') 2 . The integral over a' and over the 
ghost field C, C gives 

(det'A ^ (6.32) 

which cancels the overall factor in (6.28). 

Before proceeding, we just note that due to gauge invariance we can diagonalize ctq so that 
it takes values in the Cartan subalgebra. This introduces the usual Vandermonde factor in the 
integral over do, namely 

ct>0 
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where a denote the roots of the Lie algebra g, and a > are the positive roots. Using the Cartan 
decomposition of g, we can write as 

A lt = Y l ^X a + h li (6.34) 

a 

In this equation, X a are representatives of the root spaces of G, normalized as 

Tr(X a Xp) = 5 a+ [s, (6.35) 

where 5 a+ j3 is one if a + f3 = 0, and zero otherwise. In (6.34), is the component of along 
the Cartan subalgebra. Notice that this part of An will only contribute a co-independent factor 
to the one- loop determinant, so we will ignore it. We have 

[a ,A^] =J2<x{<To)A%X a (6.36) 

a 

and a similar equation for A. Plugging this into the action, we can now write it in terms of 
ordinary (as opposed to matrix valued) vectors and spinors 

\ j Vsd 3 * {f A 7* (" A + «(^o) 2 ) K + A - Q (V V M + ia(ao) - A") . (6.37) 

We now have to calculate the determinants of the above operators. The integration over the 
fluctuations of the gauge field is restricted, as in the Chern-Simons case, to the vector spherical 
harmonics. Using the results (4.14), (4.15), we find that the bosonic part of the determinant is: 



det(bosons) = \\ [J ((n + l) 2 + a(a ) 2 ) 2 " (n+2) . (6.38) 



>2\ 2n(n+2) 

a n=l 

For the gaugino, we can use (A. 66) to write the fermion determinant as: 

n(n+l) 



det(fermions) = \\ £J f (n + \a(ao)){—n — 1 + ia(cro)) ] , (6.39) 

a n=l ^ ' 



and the quotient gives 



ganger i _ TT TT ( n + ia(<To))"( n+1 >(-n ~ 1 + iafro))"^ 1 ) 

W ° J " y M ((n + l) 2 + a(aofy^ 

y „ = i + ia(a ))( n - 1 )( n+1 )(n + 1 - ia(a ))™(™+ 2 ) ' 



(6.40) 



up to a do-independent sign. We see there is partial cancellation between the numerator and the 
denominator, and this becomes: 

7 gaugc r i _ TT TT {n + \a{(J )) n+l TT TT (n 2 + a(a ) 2 ) w+1 

1 - 1 ° op[ ° J " A a A M (« - -(-o))- 1 " 11 11 (n* + a(a ) 2 )^ 

(6-41) 

= nn(- 2 +^o) 2 ) 2 , 

a >0n=l 



- 39 - 



where we used the fact that the roots split into positive roots a > and negative roots —a, 
a > 0. We finally obtain 



\n=l / a>0n=l 



2\ 2 



We can regularize this infinite product with the zeta function. This will lead to a finite, numerical 
result for the infinite product 

oo 

I] n\ (6.43) 

n=l 

On the other hand, we can use the well-known formula 



^M = nfl + ^) (6-44) 

n=l x ' 

to write 

ZfZ> , lWU (^m 2 , ,0.45, 

where the proportionality factor is independent of a$. We conclude that the localization of the 
vector multiplets leads to a total contribution to the partition function 

[dp [J (2smh (a (|))) 2 e"^^ (6.46) 
where we defined the convenient coupling, 

9s= 2 4 (6-47) 



and we wrote 

" 2vr' 



ao = ^, (6.48) 



where fi takes values in the Cartan subalgebra. 

6.3 Localization in Chern Simons matter theories: matter sector 

Let us now consider the matter sector. We will follow the computation in [48], which simplifies a 
little bit the original computation in [61]. As shown in (2.42), the matter Lagrangian is in itself 
a total superderivative, so we can introduce a coupling t in the form 

-tS mattcr . (6.49) 

By the by now familiar localization argument, the partition function is independent of t, as long 
as t > 0, and we can compute it for t = 1 (which is the original case) or for t — > oo. We can also 
restrict this Lagrangian to the localization locus of the gauge sector. The matter kinetic terms 
are then 

=g^d^d u <t) + fait + 2l ^~ l) 4>a Q <j) + A(2 r 7 A) #, 

A — 2 - 

C^ = - it^^d^ip + iipaoip ipip. 
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The real part of the bosonic Lagrangian is positive definite, and it is minimized (and equal to 
zero) when 

= 0. (6.50) 

Like before, in the t — > oo limit, only quadratic terms in the matter fields contribute to the 
localization computation. In particular, there is no contribution from the superpotential terms 
involving the matter multiplets, like (2.43). After using (A. 43) and (A. 61), we find that the 
operators governing the quadratic fluctuations around this fixed point are given by the operators 



= ^ {4L 2 - (A - im )(A - 2 - im )} , 

= 1 {4L • S + im + 2 - A} . 

r 

Their eigenvalues are, for the bosons, 

A (n) = r" 2 (n + 2 + ira - A)(n - ira + A), n = 0, 1, 2, • 
with multiplicity (n + l) 2 , and for the fermions 

\j,{n) = r _1 (n + 1 + ira - A), r _1 (-n + ira - A), n = l,2, 
with multiplicity n(n + 1). We finally obtain, after setting r = 1, 



IdetA 



dctA, 



n 

m>0 



(m + 1 + ircio - A) m ( m+1 )(m - ira + A) m ( m+1 ) 
(m + 1 + iroo - A) m2 (m - 1 - ira + A)™ 2 



and we conclude 



7 matterr i TT ^+ 1 ~ A + irCTp \ 



(6.51) 
(6.52) 
(6.53) 

(6.54) 
(6.55) 



m>0 



As a check, notice that, when A = 1/2 and ao = 0, we recover the quotient of determinants (4.45) 
of the free theory. The quantity (6.55) can be easily computed by using ^-function regularization 
[55, 48]. Denote 

z = 1 - A + ira (6.56) 



and 

We can regularize this quantity as 

e(z 

On the other hand, 

oo 

E 



i{z) = log Z^[a ]. 



d_ 

ds 



y , 

n ^— ' V (m + z) s (m — zV' 

u m=l 



(6.57) 



(6.58) 



m 



m 



^(m + z)' 5 (m — zY 



(h(s - 1, z) - z( H (s, z) - ( H (s - 1, -z) - z( H {s, -z), (6.59) 



where 



m=0 



(m + z) 6 



(6.60) 
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is the Hurwitz zeta function. Using standard properties of this function (see for example [80]), 
one finally finds the regularized result 

l{z) = -z log (1 - e 2 ™) + i (irz 2 + -U 2 (e 2 * iz )) - ^. (6.61) 

2 \ TT J 12 

As a check of this, notice that 




(6.62) 



in agreement with (4.44). 

There is an important property of £(z), namely when A = 1/2 (canonical dimension) one 

has 

X - (£(z) + £(z*)) = - 1 - log (2 cosh(vrr C 7 )) . (6.63) 

To prove this, we write 

z = \ + iO, (6-64) 

and we compute 

\ (£(z) + £(z*)) = -^log(2coBh(7rf)) + \tt\Q 2 + g + ±- (u 2 (-e~^ e ) + Li 2 (-e 2 ^)) . (6.65) 
After using the following property of the dilogarithm, 

Li 2 (-x) + U 2 {-x~ l ) = ~~\ (log(^)) 2 , (6.66) 

we obtain (6.63). 

When the matter is in a self-conjugate representation of the gauge group, the set of eigen- 
values of do is invariant under change of sign, therefore we can calculate the contribution of such 
a multiplet by using (6.63). We conclude that for such a matter multiplet, 

^ m £ t oTN=n( 2cosh ^) 1/2 - ( e - 67 ) 

where we set r = 1 and we used the variable /i in the Cartan defined in (6.48). The product is 
over the weights A of the representation of the matter multiplet. For general representations and 
anomalous dimensions, one has to use the more complicated result above for £{z). 

6.4 The Chern Simons matrix model 

As a first application of the results of localization, let us consider pure supersymmetric Chern- 
Simons theory, defined by the action (2.21). If we don't add matter to the theory, the fields D, 
a and A, A are auxiliary and they can be integrated out. In other words, supersymmetric Chern- 
Simons theory on S 3 should be equivalent to pure (bosonic) Chern-Simons theory. There is 
however an important difference: in super-Chern-Simons theories with at least N = 2 supersym- 
metry, there is no renormalization of the coupling k due to the extended supersymmetry [64]. 
The localization argument developed above says that the partition function of Chern-Simons 
theory on S 3 with gauge group G should be proportional to the matrix model (6.46): 

Z CS (S 3 ) oc fdn J] (2sinh^) 2 (6.68) 



- 42 - 



where we regard \i as a weight and we use the standard Cartan-Killing inner product in the 
space of weights. For example, in the case of G = U(N), if we write \x and the positive roots in 
terms of an orthonormal basis e$ of the weight lattice, 



A? 



H = ^2,Hi<%, a ij = e i -e j , i < j, (6.69) 
i=l 

we find 

/N , v 2 

nd^n( 2sinh/fi T^ e-^^-^. (6.70) 
i=l i<j ^ ' 

The proportionality constant appearing in (6.68) should be independent of the coupling constant 
k, and it is only a function of N. The matrix model (6.46) is a "deformation" of the standard 
Gaussian matrix model. It has a Gaussian weight, but instead of displaying the standard Van- 
dermonde interaction between eigenvalues (6.33) it has a "trigonometric" deformation involving 
the sinh. This interaction reduces to the standard one for small a ■ /i, which corresponds in the 
U (N) case to a small separation between eigenvalues. 

The matrix model (6.46), with a sinh kernel, was first introduced in [72]. It was later rederived 
using geometric localization techniques in [13], and abelianization techniques in [17]. As we have 
seen following [61], it can be derived in an elegant and simple way by using supersymmetric 
localization. Actually, the matrix integral appearing in the r.h.s. of (6.68) can be calculated in 
a very simple way by using Weyl's denominator formula, as pointed out in for example [4]. This 
formula reads, 

e(w)e w( - p > u = [J 2 sinh ^ • ( 6 - 71 ) 

In this formula, W is the Weyl group of G, e(w) is the signature of w, and p is the Weyl vector, 
given by the sum of the fundamental weights. Using this formula, the matrix integral reduces to 
a sum of Gaussian integrals which can be calculated immediately, and one finds 

(det(C)) 1/2 (27r 5s ) r/2 |W|e^ 2 ^ e(w)^'^ p \ (6.72) 

where C is the inverse matrix of the inner product in the space of weights (for simply connected 
G, this is the Cartan matrix), and r is the rank of G. Using again Weyl's denominator formula 
we find, 

£ eWe ^-(p) = ilA + | "Q 2sin (6 _ 73) 

where |A+| is the number of positive roots of G. The matrix integral then gives, 

(det(C)) 1 / 2 (2vrf \W\ ^^e^oy TJ 2sin (6 . 74) 

r, (I 

where we have used Freudenthal-de Vries formula 

P 2 = ^d G y. (6.75) 

The result (6.74) is indeed proportional to the partition function of Chern-Simons theory on S 3 , 
and we can use the result to fix the normalization, iV-dependent factor in the matrix integral. 
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Let us particularize for G = U(N). In this case, one has to take as C the identity matrix, and 
the value of p 2 is the same as for SU(N). Then, the matrix integral is 



3=1 



N r / ■ \ -i N—j 

IT J \ 



2 sin 



(6.76) 



k 

which is indeed proportional to (4.28) (after changing k — > k — N), up to an overall factor 

r^(2vr) Jv iV!et^ 2 -i). (6.77) 



The phase appearing here depends on k, and it has the right dependence on k, N to be understood 
as a change of framing of S 3 in the result (4.28). We can now use the above result to fix the 
normalization in the matrix model describing supersymmetric Chern-Simons theory, and we find 



i=l i<j \ ' 

We will refer to this model as the Chern-Simons matrix model. Later on we will study its large 
N limit. 

6.5 The ABJM matrix model 

Let us now consider the matrix model calculating the partition function on S 3 of ABJM theory, 
or rather its generalization [5] to the gauge group U (N\) x U (N2). The contribution of the vector 
multiplets gives in the integrand 

II (2sinh / ^)V^-^ J] (2 S inh" a ^ b ) 2 e 4^ 1 ^ 2 , (6.79) 

l<i<j<Ni ^ ' l<a<b<N 2 ^ ' 

where the opposite signs in the Gaussian exponents are due to the opposite signs in the levels. 
Since there are four hypermultiplets in the bifundamental representation, we have an extra factor 
due to (6.67), 

!=1«=1 ^ ' 

The normalization of the matrix model can be fixed by using the normalization for the Chern- 
Simons matrix model (6.78), and by comparing to the perturbative one-loop result. In this way 
we find, 

— \(n}- n?) ,JVi, JVj, , , 



1 2 J i=l a=l l<i<j<Ni V V 7 7 

II (2 smh ) 2 J] (2 cosh ) " 2 e-^fe 



l<a<b<JV 2 



(6-81) 

This model is closely related to a matrix model that computes the partition function of Chern- 
Simons theory on lens spaces L(p, 1), in particular to the model with p = 2. These models were 
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introduced in [72], and the case p = 2 was extensively studied in [4]. The matrix integral for 
p = 2 is given by, 



™>>=4^/n^nt n 

1 2 17 i=l a=l l<i<j<JVi V V 7 7 

x n ( 2 sinh (^i^) ) 2 n ( 2 c ° sh (^i^) ) 2 e ^ (Ei M?+Ea vl) ■ 

\<a<b<N 2 ^ ^ ' ' i,a ^ ^ 77 

(6.82) 

We will refer to this matrix model as the lens space matrix model. It turns out [75] that the 
partition functions (6.81) and (6.82) are related, order by order in perturbation theory, by the 
analytic continuation 

N 2 -> -N 2 . (6.83) 

We will show this explicitly in the analysis of the planar limit (indeed, we will study the planar 
limit of (6.82) rather than (6.81)). It can be shown that (6.81) is the super-matrix model version 
of (6.82) [75, 35], and this leads to the relation (6.83) between both matrix models. 

In order to derive the interpolating function for (1.1), we just have to compute the planar free 
energy of the matrix model (6.81) for N± = N 2 = N, and for any value of the 't Hooft coupling 
A. The calculation of exact planar free energies of matrix models is a classical problem which 
was first solved in [18] for a simple class of matrix models. We will now review the standard 
techniques to do that, which we will then generalize to the matrix models appearing in ABJM 
theories. 

7. Matrix models at large N 

In this section we will focus on conventional matrix models, and in the next section we will 
use the same techniques (and the same formulae) to analyze the matrix models appearing in 
supersymmetric Chern-Simons-matter theories. A more detailed treatment of matrix models in 
the large expansion, as well as a complete list of references, can be found in [29, 74, 39]. 

7.1 Saddle-point equations and one-cut solution 

Let us consider the matrix model partition function 

1 i r N 

* = M^/n«<A*(A)e-A^.> (7.1) 

1=1 

V(A), called the potential of the matrix model, will be taken to be a polynomial 

VW = \x 2 + E 9 i; XP ( 7 - 2 ) 

where the g p are coupling constants of the model. In (7.1), 

A 2 (A) = - Xj) 2 (7.3) 

i<j 
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is the Vandermonde determinant (6.33) for the group U(N). The integral (7.1) is typically 
obtained as a reduction to eigenvalues of integrals over the space of N x N Hermitian matrices, 
see [29, 74] for more details. We want to study Z in the so-called 't Hooft limit, in which 

g s -»• 0, N -»• oo, (7.4) 

but the 't Hooft parameter of the matrix model 

t = g s N (7.5) 

is fixed. In particular, we want to study the leading asymptotic behavior of the free energy 

F = logZ (7.6) 

in this limit. Let us write the partition function (7.1) as follows: 

N 
i=l 

where the effective action is given by 

s cS W = ~ E v(\i) + |J E lo s l A * - ( 7 - 8 ) 

i=l i<j 

We can now regard gl as a sort of h, in such a way that, as g s — > with i fixed, the integral (7.7) 
will be dominated by a saddle-point configuration that extremizes the effective action. Notice 
that, since a sum over N eigenvalues is roughly of order N, in the 't Hooft limit the effective 
action is of order 0(1), and the free energy scales as 

F(g s ,t)*g- 2 F (t). (7.9) 

Fo(t) is called the genus zero, or planar, free energy of the matrix model, and it is obtained by 
evaluating the effective action at the saddle point. This dominant contribution is just the first 
term in an asymptotic expansion around g s = 0, 

oo 

F = Y J m)g 2 r 2 - (7.io) 

9=0 

In order to obtain the saddle-point equation, we just vary S e g(\) w.r.t. the eigenvalue Aj. We 
obtain the equation 

This equation can be given a simple interpretation: we can regard the eigenvalues as coordinates 
of a system of N classical particles moving on the real line. (7.11) says that these particles are 
subject to an effective potential 

y cff (A,) = V(Xi) - | J] log |A< - Xj\ (7.12) 



-46- 



which involves a logarithmic Coulomb repulsion between eigenvalues. For small 't Hooft param- 
eter, the potential term dominates over the Coulomb repulsion, and the particles tend to be at 
a critical point of the potential: V'(x*) = 0. As t grows, the logarithmic Coulomb interaction 
will force the eigenvalues to repel each other and to spread out away from the critical point. 

To encode this information about the equlibrium distribution of the particles, it is convenient 
to define an eigenvalue distribution (for finite N) as 

1 N 

pCAH^JXA-A,)), (7.13) 

i=l 

where the Aj solve (7.11) in the saddle-point approximation. In the large N limit, it is reasonable 
to expect that this distribution becomes a continuous distribution po(X). As we will see in a 
moment, this distribution has a compact support. The simplest case occurs when po(X) vanishes 
outside a connected interval C = [a,b]. This is the so-called one-cut solution. Based on the 
considerations above, we expect C to be centered around a critical point of the potential. In 
particular, as t — > 0, the interval C should collapse to the point x*. 

We can now write the saddle-point equation in terms of continuum quantities, by using the 

rule 



1 N r 

*f£/M-> / /(A)A)(A)dA. (7.14) 



1=1 

Notice that the distribution of eigenvalues po(A) satisfies the normalization condition 



Po(X)dX = 1. (7.15) 

c 



The equation (7.11) then becomes 



where P denotes the principal value of the integral. The above equation is an integral equation 
that allows one in principle to compute po(X), given the potential V(X), as a function of the 't 
Hooft parameter t and the coupling constants. Once po(A) is known, one can easily compute 
Fo(t), since the effective action in the continuum limit is a functional of po : 

S cS (p ) = -t f dX Po (X)V(X) + t 2 [ dAdA / po(A)p (A / )log|A-A / |. (7.17) 

JC JCxC 
The planar free energy is given by 

F (t) = Seff(po)- (7.18) 

We can obtain (7.11) directly in the continuum formulation by computing the extremum of 
the functional 

S(p ,T) = S eS (p ) + r(t ^dXpo(X) - ?J (7.19) 

with respect to po- Here, T is a Lagrange multiplier that imposes the normalization condition of 
the density of eigenvalues (times t) . This leads to 

V(A) = 2t J dX' po(X') log |A - A' | + T, (7.20) 
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which can be also obtained by integrating (7.16) with respect to A. It is convenient to introduce 
the effective potential on an eigenvalue as 



V cG (X) = V(X)-2t J dA / p (A / )log|A-A / |. (7.21) 

This is of course the continuum counterpart of (7.12). In terms of this quantity, the saddle-point 
equation (7.20) says that the effective potential is constant on the interval C: 

V cS (X) = T, A g C. (7.22) 

The Lagrange multiplier V appears in this way as an integration constant that only depends on 
t and the coupling constants. As in any other Lagrange minimization problem, the multiplier 
is obtained by taking minus the derivative of the target function w.r.t. the constraint, which in 
this case is t. We then find the very useful equation 

d t F {t) = -T = -V eS (b). (7.23) 

where b is the endpoint of the cut C. 

The density of eigenvalues is obtained as a solution to the saddle-point equation (7.16). This 
equation is a singular integral equation which has been studied in detail in other contexts of 
physics (see, for example, [79]). The way to solve it is to introduce an auxiliary function called 
the resolvent. The resolvent is defined, at finite N, as 



"(p) = ii?(E^r>. ( 7 - 24 ) 




and we will denote its large N limit by coo(p), which is also called the genus zero resolvent. This 
can be written in terms of the eigenvalue density as 

u*(p) = J dX^. (7.25) 

The genus zero resolvent (7.25) has three important properties. First of all, due to the normal- 
ization property of the eigenvalue distribution (7.15), it has the asymptotic behavior 

w o(p) ~ -j P — > °°- (7.26) 

Second, as a function of p it is an analytic function on the whole complex plane except on the 
interval C, where it has a discontinuity as one crosses the interval C. This discontinuity can be 
computed by standard contour deformations. We have 

MP + V = ! dA^^- = / dA^l = P / dA^l + / dA^l, (7.27) 
Jr p + ie- x y R _ ie p — X J p - X J Ce p-X 

where C t is a contour around A = p in the lower half plane, and oriented counterclockwise. The 
last integral can be evaluated as a residue, and we finally obtain, 

co (p + ie) = P J dA^ - nipo(p). (7.28) 
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Similarly 

MP ~ k) = / = P / dA^J + nipoip). (7.29) 
One then finds the key equation 

PoW = -^7(wo(A + ie)-wo(A-ie)). (7.30) 

From these equations we deduce that, if the resolvent at genus zero is known, the planar 

eigenvalue distribution follows from (7.30), and one can compute the planar free energy. On the 
other hand, by using again (7.27) and (7.29) we can compute 



lo (p + ie) + u (p - ie) = 2P J dA^^ (7.31) 



and we then find 

co {p + ie) + cu (p - ie) = -V(p), p € C, (7.32) 

which determines the resolvent in terms of the potential. In this way we have reduced the original 
problem of computing Fo(t) to the Riemann-Hilbert problem of computing u>o(A). In order to 
solve (7.32), we write it as a sum of an analytic or regular part w r (p), and a singular part u> s (p), 

u (p) = uj r (p) + uj s {p), (7.33) 

where 

wr(p) = jV'(p). (7.34) 
It follows that the singular part satisfies 

uj s {p + \e) +uj s (p-ie) = 0, pGC. (7.35) 
This is automatically satisfied if oo s {p) has a square-root branch cut across C, and we find 

"s(p) = ~M(p)y/(p - o)(p - b), (7.36) 

where a,b are the endpoints of C, and M(p) is a polynomial, which is fully determined by the 
asymptotic condition (7.26). There is in fact a closed expression for the planar resolvent in terms 
of a contour integral [77] which reads 

„ oM 1 I^mfp^zR) 1 ' 2 , (7 . 37 ) 
yFJ 2t J c 2mp- z\(z-a)(z-b) J ' v 1 

where C denotes now a contour encircling the interval. The r.h.s. of (7.37) behaves like c + d/p + 
0(l/p 2 ). Requiring the asymptotic behavior (7.26) imposes c = and d = 1, and this leads to 

V (7.38) 
dz zV'(z) _ 2t 
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Figure 3: The contour C encircling the support of the density of eigenvalues, which can be regarded as 
a contour in the spectral curve y(p). 



These equations are enough to determine the endpoints of the cuts, a and b, as functions of the 
't Hooft coupling t and the coupling constants of the model. Equivalently, after deforming the 
contour in (7.37) to infinity, we pick a pole at z = p, which gives the regular piece, and we find 
the equation 

"*(P) = Yt V ' ip) ~ Y t M(p)V(p-a)(p-b), (7.39) 

M( P ) = } ^ V ' {Z) , 1 = . (7.40) 
1 2iriz-p^( z - a )(z-b) V ' 



where 



A useful way to encode the solution to the matrix model is to define the spectral curve of 
the matrix model by 

yip) = V'{p) - 2tu (p) = M(p)y/(p-a)(p-b). (7.41) 
Notice that, up to a constant, 



/ 



A 

dpy(p) = V cS (X). (7.42) 



If we regard uo(p)dp as a differential on the spectral curve, the 't Hooft parameter can be written 
as a contour integral 

t = j ^2fwo(p). (7.43) 
This contour on the spectral curve (regarded as a complex curve) is represented in Fig. 3. 

Example 7.1. The Gaussian matrix model. Let us now apply these results to the simplest case, 
the Gaussian model with V{z) = z 2 /2. We first look for the position of the endpoints from 
(7.38). After deforming the contour to infinity and changing z — > 1/z, the first equation in (7.38) 
becomes 

(k 1 ' 0, (7.44) 



o 2vri z 2 _ az )(l - bz) 

where the contour is now around z = 0. Therefore a + b = 0, in accord with the symmetry of 
the potential. Taking this into account, the second equation becomes: 

— -1 , 1 ==2t, (7.45) 
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and gives 

a = 2\Tt. (7.46) 

We see that the interval C = [—a, a] = [-2y/i,2y/i\ opens as the 't Hooft parameter grows up, and 
as t — > it collapses to the minimum of the potential at the origin, as expected. We immediately 
find from (7.39) 

Mp) = Y t {p~ \/p 2 -4t), (7.47) 
and from the discontinuity equation we derive the density of eigenvalues 

p (A) = t^V^-X 2 . (7.48) 

The graph of this function is a semicircle of radius 2\/t, and the above eigenvalue distribution is 
the famous Wigner-Dyson semicircle law. Notice also that the equation (7.41) is in this case 

y 2 =p 2 - At. (7.49) 

This is the equation for a curve of genus zero, which resolves the singularity y 2 = p 2 . We then 
see that the opening of the cut as we turn on the 't Hooft parameter can be interpreted as a 
deformation of a geometric singularity. 

7.2 Multi cut solutions 

So far we have considered the so-called one-cut solution to the one-matrix model. This is not, 
however, the most general solution, and we will now consider the so-called multi-cut solution, in 
the saddle-point approximation. Recall from our previous discussion that the cut appearing in 
the one-matrix model was centered around a critical point of the potential. If the potential has 
many critical points, one can have a saddle-point solution with various cuts, centered around 
different critical points. The most general solution has then n cuts (where n is lower or equal 
than the number of critical points), and the support of the eigenvalue distribution is a disjoint 
union of n intervals 

n 

C = [jC t . (7.50) 
i=i 

The total number of eigenvalues N splits into n integers N{, 

N = N x + • • • + N n , (7.51) 

where iV, is the number of eigenvalues in the interval Cj. We introduce the filling fractions 

N- I' 

tl = ~W = / dA/9 °( A )' * = !,■•• ,n. (7.52) 

Notice that 

n 

5> = i. ( 7 - 53 ) 

A closely related set of variables are the partial 't Hooft parameters 

ti = tei = g s Ni, i = l,---,n. (7.54) 
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Notice that there are only g = n — 1 independent filling fractions, but the partial 't Hooft 
parameters are all independent. 

The multi-cut solution is just a more general solution of the saddle-point equations that we 
derived above. It can be found by extremizing the functional (7.17) with the condition that the 
partial 't Hooft parameters are fixed, 

S( Po , e 1 ) = S eS (p ) + it T *( t J c dX ~ **) ' ( 7 - 55 ) 

where Tj are Lagrange multipliers. If we take the variation w.r.t. the density po(\) we find the 
equation 

V(A) = 2t^dA'/?o(A / )log|A- A'| + 1^, A e d (7.56) 
which can be rewritten as 

Kff(A) =Fi, AG Cj. (7.57) 

The planar resolvent still solves (7.32), and the way to implement the multi-cut solution is to 
require wo(p) to have 2n branch points. Therefore we have 



2n 



v Hip-Xk), (7.58) 
\ k=i 



which can be solved in a compact way by 



(2n \ ^"/^ 

In order to satisfy the asymptotics (7.26) the following conditions must hold: 

1 I dz z e V'(z) 



'J£n 



2tIc2ni U l li{z _ Xk) - 2 



£ = 0,l,---,n. (7.60) 



In contrast to the one-cut case, these are only n+1 conditions for the 2n variables Xk representing 
the endpoints of the cut. The remaining n—1 conditions are obtained by fixing the values of the 
filling fractions through (7.52) (or, equivalently, by fixing the partial 't Hooft parameters). The 
multipliers in (7.55) are obtained, as before, by taking derivatives w.r.t. the constraints, and we 
find the equation 

-du'du^- rm " r *' (7 - 61) 

which generalizes (7.23) to the multi-cut situation. 

We can write the multi-cut solution in a very elegant way by using contour integrals. First, 
the partial 't Hooft parameters are given by 

1 



, , _ 2tiv (p)dp. (7.62) 

47T1 If. 

We now introduce dual cycles T>i cycles, i = 1, ■ ■ ■ , n — 1, going from the Cj+i cycle to the Ci 
cycle counterclockwise, see Fig. 4. In terms of these, we can write (7.61) as 

dF dF 



dU dt i+ i 2 



X - I 2tu (p)dp. (7.63) 

JT>i 
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Figure 4: A two-cut spectral curve, showing two contours C\,2 around the cuts where Nx,2 eigenvalues 
sit. The "dual" cycle T>\ goes from C 2 to C\. 



8. Chern— Simons and ABJM matrix models at large N 

In this section we will solve the matrix models (6.78) and (6.81), in the planar limit, by using 
the saddle-point techniques of the last section. 

8.1 Solving the Chern Simons matrix model 

Let us first consider the Chern-Simons matrix model (6.78). The analogue of the saddle-point 
equation (7.11) is 



1 



Hi - 



The form of the r.h.s. suggest to define the resolvent as 

JV 

z- Hi 

com i 

The large ./V limit of (8.1) gives then 



{z) =9s(y2 coth \~2^) 



z = \ ( u o( z + ie ) + ^0(2 - ie)) • 
The planar resolvent satisfies the boundary conditions 

ojq(z) ~ ±t, z — > ±00, 

where 



13) 



J.4) 



t = g s N (8.5) 
is the appropriate 't Hooft parameter for this model. Let us define the exponentiated variable 

Z = e z . (8.6) 

In terms of the Z variable, the resolvent is given by 



N 



dZ I ^ dZ 

(z)dz = -t— + 2g s (}2 



,i=l 



Z - e^ 



(8.7) 



Prom this resolvent it is possible to obtain the density of eigenvalues at the cuts. In the planar 
approximation, we have that 



uj (z) = -t + 2t 



Z-^ 



dfi, 
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where po(fi) is the density of eigenvalues, normalized in the standard way 

po{n)dn = 1. 



(8.9) 



The discontinuity argument which gave us (7.30) in the standard matrix model case tells us now 
that 

p {X)dX = --±-^(uo{X + ie)-GJ (X-k)), XGC. (8.10) 

47Tlt A 

Let us solve explicitly for ujq by using analyticity arguments, following [47]. We first construct 
the function 

g{Z) = e ^ 2 + Ze-" o/2 . (8.11) 
This function is regular everywhere on the complex Z plane. Indeed, we have 

g(Z + ie) = e^ z +*)/ 2 + Ze ~^ z+it ^ 2 = Zq~^ z -' 1 ^I 2 + e~^ z ~ ie ^ 2 = g{Z - ie), (8.12) 

so it has no branch cut. The boundary conditions for this function, inherited from (8.4), are 

lim g(Z) = e~ t/2 Z, lim g(Z) = e~ t/2 . (8.13) 

These conditions are solved by 



g(Z) = e- t / 2 (Z + l), 
and we can now regard (8.11) as a quadratic equation that determines ojq: 



\ (g(Z) - ^g*{Z)-Az) 



uj q {Z) = 2\og 

From this resolvent we can determine immediately the density of eigenvalues, 

1 



(8.14) 



(8.15) 



Po( x ) 



irt 



tan 



cosh 



2 /x' 



cosh (|) 



3.16) 



supported on the interval [—A, A] with 

,4 = 2 cosh" 1 (e*/ 2 ) . (8.17) 

The result (8.15) can be also obtained directly from (7.37), see [74]. 
8.2 Solving the ABJM matrix model 

Let us now solve the matrix model we are interested in, namely (6.81). The saddle point equations 
for the eigenvalues /ij, u a are 

Ni N 2 



— = coth : , 

9s ^ 2 ^ 2 

N 2 Ni 

Va = ^ CQth v a - I'b \- , ,,,, "a - Ih 

9 S l,J-r. 



a=l 

tanh 



(8.18) 



i=l 
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We will solve instead the saddle-point equations 

. m N 2 
* =^ 2. coth ~ + ^ 2. tanh ^— , 

i^i a=l 

(8.19) 

. N 2 , Ni 

= ^ 2. coth — g— + ^ !L tanh — g— • 

b^a i=l 

where 

ti = 5siVi (8.20) 

are the partial 't Hooft parameters for this model. Clearly, from the solution to (8.19) we can 
recover the solution to (8.18) by simply performing the analytic continuation 

t 2 -»• -t 2 (8.21) 

in the solution. The equations (8.19) are the saddle-point equations for the lens space matrix 
model (6.82), and the analytic continuation (8.21) is just the planar version of the relation (6.83) 
mentioned before. 

The total resolvent of the matrix model, oj(z), is defined as [47] 

u,(z) = 9s coth ^ + 9s (f: tanh ^ . (8.22) 



In terms of the Z variable, it is given by 

uMd, = - t f + 2 3 , (t ™ ) + 2 9 , <£ ^ ) ■ (8.23) 



Z -eK / \ ^ Z + e"» 

i i=l / \a=l 



where t = ii + ^2, and it has the following expansion as Z — > oo 

iVi jV 2 \ 

z -l. 



\i=l a=l 



From the total resolvent it is possible to obtain the density of eigenvalues at the cuts. In the 
planar approximation, we have that 

wo(z) = -t + 2ti / Pl (fi)-^—d^ + 2t 2 [ p 2 {v)^—dv, (8.25) 
JCi z - eM ic 2 z + e 

where pi(/i), P2(^) are the large iV densities of eigenvalues on the cuts Ci, C2, respectively, 
normalized as 

/ pi(M)d/i= / p 2 {v)dv = l. (8.26) 

The standard discontinuity argument tells us that 

1 dX 

Pl (X)dX = --r—^- (uq(X + ie) - oj (X - ie)) , X € d, 
47rici A 

(8.27) 

p 2 (y)dy = (u (Y + ie) - MY ~ ie)) , ^^2. 
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Let us now find an explicit expression for the resolvent, following [47]. First, notice that it 
can be split in two pieces, 

u(z) = u { - 1) (z)+J 2) {z + m), (8.28) 

where 

/ Ni , 



u>W(z) = 5s ( ^coth 
\i=i 



N2 

s \ ^2 coth 



1.29) 



are just the resolvents of the Chern-Simons matrix model (8.2). In fact, it is easy to see that 
the matrix model (6.82) is equivalent to a Chern-Simons matrix model for N variables U{, 
£ = !,••• ,N, where N\ variables 



are expanded around the point z = 0, and N2 = N — N\ variables 



^JVi+a = ivr + i/ a , a = 1, • 



(8.30) 



5.311 



are expanded around the point z = ivr. At large iVi 2 it is natural to assume that the first set of 
eigenvalues will condense in a cut around z = 0, and the second set will condense in a cut around 
z = 7ri. It follows that u^ l \z) will have a discontinuity on an interval [— A, A], while will 
have a discontinuity on an interval [—B,B]. When g s is real, these cuts occur in the real axis, 
and the two cuts in the total resolvent are separated by iir (see Fig. 5). 



V 



i<5 



m-B 



C, 



L£ 



7ri + B 
Ci 



A 



C, 



■1/6 



2? 



Ci 



1/a a 



Figure 5: Cuts in the z-plane and in the Z-planc. 
The saddle-point equations (8.19) become then, at large N, 



z =- (u)q(z + ie) + u (z - ie)) , z £ [-A, A], 

z =- (uq(z + ivr + ie) + u>o(z + ivr - ie)) , z G [-B, B]. 
It follows that the function 

f{Z) = e* (e^° + ZV^ ) 
is regular everywhere on the complex plane and has limiting behavior 



lim f(Z) = Z< 

Z—>oo 



lim f(Z) = l. 



(8.32) 

(8.33) 
(8.34) 
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The unique solution satisfying these conditions is 



f(Z) = Z 2 -(Z + l, (8.35) 

where ( is a parameter to be determined. Solving now (8.33) as a quadratic equation for e w ° 
yields, 

uo(Z) = log^ [f(Z) - y/p(Z)-4e*Z*\ ) . (8.36) 

Notice that e^° has a square root branch cut involving the function 

a(Z) = f(Z) - 4e 2t Z 2 = (Z - a) (Z - l/a) {Z + b) (Z + 1/b) (8.37) 

where a 1 * 11 , — 6 ±1 are the endpoints of the cuts in the Z = e z plane, see Fig. 5. We deduce that 
the parameter ( is related to the positions of the endpoints of the cuts by the relation 

^^(H- 6 -!)- (8 - 38) 

and we also find the constraint 

i( o+ i +6+ £) =e ' (s - 39) 

Once the resolvent is known, we can obtain both the 't Hooft parameters and the derivative 
of the genus zero free energy in terms of period integrals. The 't Hooft parameters are given by 

U = — I u> (z)dz, i = l,2. (8.40) 
4vri Jc- 

This is the same equation that (7.62), since the resolvent (8.22) has an extra factor 2t as compared 
to (7.24). It is easy to see, by using the same techniques that we used in section 7, that the 
planar free energy Fq satisfies the analogue of (7.63) after taking into account this extra factor, 



ds 

where 



- vrit = -\ (f oj (z)dz, (8-41) 

s = \{ti-t 2 ) (8.42) 

and the T> cycle encloses, in the Z plane, the interval between —1/6 and l/a (see Fig. 5). The 
extra term — -zrit in (8.41) is due to the fact that ujq{z)&z has a nonzero residue at z = 0. 

The above period integrals are hard to compute, but their derivatives can be easily found 
by adpating a trick from [19]. One obtains, 

= <P , ^==±^ T^K(k), (8.43 

d( 4vri J ^(Z 2 - (Z + l) 2 - 4e 2 *Z 2 vr(l + ab) v h v 1 

Cl,2 

where K(k) is the complete elliptic integral of the first kind, and its modulus is given by 

(« 2 - 1)(& 2 - 1) _ , (« + b\ 2 (8 . 44) 



(l + a&) 2 \l + ab 
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Likewise for the period integral in (8.41) we find 



where the complementary modulus k' = y/l — k 2 is given by 

k> = ^±. (8.46) 
1 + ab y ' 

8.3 ABJM theory and exact interpolating functions 

We will now analyze the exact planar results for the ABJM matrix model in the simplest case, 
namely the one corresponding to the original ABJM theory with gauge groups of the same rank 
N. If we recall the definition of g s (6.47) in terms of the level k, we find 

N , 
t x = -t 2 = 2vri— = 27riA, (8.47) 

where A is the 't Hooft coupling of ABJM theory defined in (1.2). If we think about ti % as "mod- 
uli" parametrizing the space of complex 't Hooft couplings, the ABJM theory of [6] corresponds 
to a real, one-dimensional submanifold in this moduli space. We will call this submanifold the 
ABJM slice. In this slice, the total 't Hooft parameter of the matrix model vanishes: t = 0. 
The theory has only one parameter, A, which should be related to the only parameter (8.38) 
appearing in the resolvent. It follows from (8.38) and (8.39) that 

a+- = 2 + C, b+]- = 2-C, (8.48) 
a b 

The derivative (8.43) can be expressed in a simpler way by using appropriate transformations of 
the elliptic integral K(k). Let us consider the elliptic moduli 

l — k' , ki , 

= k2 = 'W i8A9 > 

One has that (see for example [45], 8.126 and 8.128) 

K(k) = (1 + fci) J KT(fc 1 ) = i±^tf(*2), (8.50) 

and we can write 



vr(l + a6) w vr(l + a6) k[ y ' ir V (a + 6)(1 + ab) 
Notice that 

kl ~ (a + l)(6+l)' k2 ~ A(a + b)(l + ab) ^ 
In the ABJM slice we have 



, 2 C 2 / ab 1 



fe2 "l6' V (a + b)(l + ab)-2 (8 - 53) 
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and 



dA 



4tt 2 



K 



(8.54) 



Since K(k) depends only on k 2 , t follows from this equation that, if we want A to be real (as it 
should be in the ABJM theory), ( has to be pure imaginary, and we can write 



so that (8.54) becomes 



C = IK, 

dA 1 



K G 



dn 47T 2 ' 

This can be integrated explicitly in terms of a hypergeometric function [75] 

^2- 



A(«) 



8vr 



3-^2 



111 3 

2'2'2' '2' 



ft 

16 



(8.55) 
(8.56) 

(8.57) 



where we have used that A = when k = (in this limit, the cut [a, 1/a] collapses to zero 
size, and the period t\ vanishes). (8.57) gives the relation between A and ( (or k). This closed 
expression makes it possible to perform an analytic continuation for k 1, where A behaves as 



A(«) 



log 2 (*Q 
2tt 2 



1 



This suggests to define the shifted coupling 

A = A 



1 

24' 



(8.58) 



(8.59) 



Notice from (5.80) that this shift is precisely the one needed in order for A to be identified with 
Q/k at leading order in the string coupling constant. The relationship (8.58) is immediately 
inverted to 

1 



,7iV2A 



l + O 



,e 



-2ttV2A 



(8.60) 



Let us now consider the genus zero free energy. Its second derivative w.r.t. s, evaluated at 
t = 0, can be calculated as 



d 2 F n 



da 



dX 




t=o d( 


t=o ( 



'dti 



(8.61) 



Like before, we will use the transformation properties of the elliptic integral K(k') to write (8.45) 
in a more convenient way. Prom (8.49) we deduce 



k[ 



2Vk> 



1 + k r 

and we have, using again [45], 8.126 and 8.128, 



1 



k[ 



K(k') 

In the ABJM slice we find, 



1 



^(K{k' 2 ) + iK(k 2 )). 



81 



t=0 



K' 



+ '\K 



(8.62) 



(8.63) 



(8.64) 
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so that 



d 2 F 



ds 2 

We conclude that, in the ABJM theory, where s = 27riA, 



K' {—) 

= -it )Af - 7T1. (8.65) 

t=o K \T) 



^(^=4^-^11 + 4^. (8.66) 
A further integration leads to the following expression in terms of a Meijer function 

^)^^(H;A|4)^3^(^.^.|4)- 

This is, indeed, the exact interpolating function we were looking for! To see this, we can expand 
it at weak coupling as follows: 

d x F (X) = -8vr 2 A (log - l) + + O (A 5 ) . (8.68) 

After integrating w.r.t. A and multiplying by 

9s 2 = ~ 2 , (8-69) 

we find that the first term exactly reproduces the weak-coupling answer (4.49). The comparison 
with the weak coupling expansion also fixes the integration constant, 

F (A) = f dX'dyF (X'). (8.70) 
J 

To study the strong-coupling behavior, we can now analytically continue the r.h.s. of (8.67) to 
k = 00, and we obtain 

9\Fq{X) = 2-7T 2 log k + —y- 4^3 fl,l,^;2,2,2;-^ . (8.71) 



K 2 ^'^'2' 2' ' ' ' k\ 
After integrating w.r.t. A and using the shifted coupling A defined in (8.59), we find, 

F (A) = !^A 3 / 2 + Y.^^h ( -±=) (8-72) 

where fe(x) is a polynomial in x of degree 2£ — 3 (for i > 2). If we multiply by g~ 2 , we find 
that the leading term agrees precisely with the prediction from the AdS dual in (1.1). The series 
of exponentially small corrections in (8.72) were interpreted in [32] as coming from worldsheet 
instantons of type IIA theory wrapping the CP 1 cycle in CP 3 . This is a novel type of correction 
in AdS4 dualities which is not present in AdSs spaces, see [21] for a preliminary investigation of 
these effects. 



9. Further developments and results 

In these lectures our goal has been to find the interpolating function for (1.1) in ABJM theory, 
but many other results can be obtained with the ideas and techniques reviewed here. In this last 
section we will give a brief summary of more recent results along these lines. 
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Figure 6: Comparison of the exact result for d\Fo(A) given in (8.67), plotted as a solid blue line, and 
the weakly coupled and strongly coupled results. In the figure on the left, the red dashed line is the 
supergravity result given by the first term in (8.72), while in the figure on the right, the black dashed line 
is the Gaussian result given by the first two terms in (8.68). 

9.1 Wilson loops and higher genus corrections 

One of the main results of the localization techniques developed in [82], in the context of four- 
dimensional gauge theories, was a proof of the conjecture of [38, 31]. According to this conjecture, 
the vacuum expectation value of 1/2 BPS Wilson loops can be calculated by a matrix model. 
Similarly, the original motivation of [61] was to provide a matrix model calculation of Wilson 
loops in ABJM theory, generalizing [82]. One can indeed show that both the vev of the 1/6 BPS 
Wilson loops of [34, 86, 23], as well as of the vev of the 1/2 BPS Wilson loop of [35], can be 
computed as correlators in the matrix model (6.81). Using the resolvents that we derived in these 
lectures, one obtains exact interpolating functions for the vevs of these Wilson loops [75, 32], 
which at strong coupling are in agreement with the AdS predictions. 

When Ni = N2, i.e. for the original ABJM theory, the matrix model (6.81) can be explicitly 
solved to all orders in the 1/N expansion [32] by using techniques from topological string theory 
(see [46] and references therein). This is a remarkable fact, since it amounts to solving for the 
free energy of the dual type IIA superstring theory at all genera, and it makes it possible to 
address interesting nonperturbative questions [33]. 

9.2 More general Chern Simons matter theories 

In these lectures we have only considered the ABJM theory, but there are many other super- 
symmetric Chern-Simons-matter theories with proposed large N AdS duals. According to the 
AdS4/CFT3 correspondence, the free energy on the three-sphere of these theories should behave, 
in the planar limit and strong coupling, as (5.89). This probes the volume of the compactification 
manifold in M-theory, which can be a non-trivial function of other parameters of the model, and 
suggests a strategy to make precision tests of the AdS/CFT correspondence: compute the planar 
free energies of the matrix models associated to more general Chern-Simons-matter theories, 
and see if they display the behavior (5.89) at strong 't Hooft coupling. 

However, these matrix models become increasingly harder to solve as we move away from 
the ABJM theory, even in the planar limit, and exact results for the resolvent are only known for 
the Gaiotto-Tomasiello theories [92] and for theories with matter in the fundamental [26]. Even 
when the exact resolvent is known, extracting the strong coupling behavior might be difficult. 
For this reason, it is important to have techniques that make possible to find the strong coupling 
limit of the free energy, even without having an exact expression. In [26] a geometric technique 
based on ideas of tropical geometry was proposed to derive in an easy way the strong coupling 
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behavior from the resolvent, and the behavior (5.89) was tested for a one-parameter family of 
tri-Sasaki-Einstein manifolds Xj dual to N = 3 theories with fundamental matter. Moreover, in 
[52], a very powerful method was introduced which makes possible to analyze the matrix models 
describing Chern-Simons-matter theories directly in the strong coupling limit, without the need 
to solve first for the resolvent. This allowed to test (5.89) for a large class of N = 3 theories. 

9.3 Extremization of Z 

Perhaps the most interesting spinoff of these developments, from the QFT point of view, is the 
idea that the free energy of QFTs on the three-sphere should play the role, in three dimensions, 
of Zamolodchikov's c-function in two dimensions, and of the anomaly a-coefficient in four di- 
mensions: it has been conjectured that this quantity decreases along RG trajectories and it is 
stationary at RG fixed points. This has been dubbed the "F-theorem" in [56] . The first evidence 
for this result comes from the extension of the localization computation of [61] to theories with 
N = 2 supersymmetry [55, 48]. In this case, the anomalous dimensions of the matter fields A 
are not canonical, and the partition function is itself a function of them. It was shown in [55] 
that the partition function was extremized as a function of these anomalous dimensions, giving 
an efficient way to calculate them. This result is similar to the four-dimensional result that the 
R-charges extremize the anomaly a-coefficient [54], and it leads naturally to the "F-theorem" 
conjecture. The extremization property proved in [55] has been applied and tested in various 
situations [76, 24, 56, 9], although the validity of the F-theorem is still conjectural and possible 
counterexamples have been discussed in [56, 81]. 

9.4 Other extensions and applications 

The localization techniques of [61] have been used as well to test field theory dualities in three- 
dimensional supersymmetric gauge theories [62, 63, 57, 60, 96], by showing that their partition 
functions agree as functions of the various parameters. They have also been used to compute 
supersymmetric indices in § 2 x S 1 [65, 53, 69] and to study Chern-Simons-matter theories in 
other three-manifolds [49, 43, 59]. 

In conclusion, the localization ideas of [82, 61] give us new ways of computing exact quantities 
in gauge theories where matrix model techniques play a crucial role, and we anticipate many 
interesting developments coming from this line of research. 
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A. Harmonic analysis on S 3 
A.l Maurer Cartan forms 

We will first introduce some results and conventions for the Lie algebra and the Maurer-Cartan 
forms. The basis of a Lie algebra g satisfies 

[T a ,T b ] = f abc T c . (A.l) 

If g G G is a generic element of G, one defines the Maurer-Cartan forms uj a through the equation 

g- 1 dg = Y,TaOJ a , (A.2) 

a 

and they satisfy 

duJa + ^fabc^b A w c = 0. (A.3) 

This is due to the identity 

d(g- 1 dg)+g- 1 dgAg- 1 dg = 0. (A.4) 
Let us now specialize to SU(2). A basis for the Lie algebra is given by: 



~ 2 



(A.5) 



Explicitly 

^Co)< (A6) 

The structure constants are 

fabc = -e a bc- (A. 7) 

Any element of SU(2) can be written in the form 

We parametrize this element as (see for example [93]) 

II tl \a\ • a h + t 3 t 2 -t 3 + 7T 

M = cos— , \/3\=sm— , Arga = — - — , Arg/?= , (A.9) 

where U are the Euler angles and span the range 

< h < ir, < t 2 < 2vr, -2vr < t 3 < 2ir. (A.10) 

The general element of SU(2) will then be given by 

/ cos(ti/2)e i (* 2 +* 3 )/ 2 isin(t 1 /2)e i (*2-*3)/2\ 
g - u(h,t 2 .t 3 ) - ^ isin(£l/2)e i(-t 2 +* 3 )/2 CO s(t 1 /2)e- i ^+*3)/2j (A n) 

= «(*2,0,0H0,ti,0)«(0,0,t 3 ). 

We then have 



_ l _\_( dt 3 + costidt 2 e lt3 (dti + idi 2 sin t{j \ , , 

5 5 ~ 2 Ve^dii - idt 2 sin h) -dt 3 - cos ^dta ' ' ' 



- 63 - 



Therefore, 

oj\ = cos tz&ti + sin ts sin tidi 2 , 

UJ2 = sin t^dti — cos £3 sin iidt 2 , (A. 13) 

W3 = COS tidf 2 + d^3, 

and one checks explicitly 

duJ a = -j^abcUb A UJ C , (A.14) 

as it should be according to (A. 3). 
A. 2 Metric and spin connection 

The metric on SU (2) = § 3 is induced from the metric on C 2 

ds 2 = r 2 ^d|a| 2 + |a| 2 dArga 2 + d|/3| 2 + |/?| 2 dArg/3 2 ^ , (A.15) 
where r is the radius of the three-sphere. A simple calculation leads to 



ds 2 = T — ( dt\ + dtj + dt§ + 2 cos ti dt 2 dt 3 ) , (A. 16) 



with inverse metric 

4 A 

G' 1 = -2 esc 2 ti -cotticscii) (A.17) 

r \0 — cot ti CSC tl CSC 2 ti 

and volume element 

(detG)V 2 = r3sin l \ (A.18) 
8 

The volume of S 3 is then 

f (det G) 1/2 dti dt 2 dt 3 = 2vr 2 r 3 (A. 19) 

JSU{2) 

which is the standard result. The only nonzero Christoffel symbols of this metric are 

T\ z = ^smt u r 2 3 = r 3 2 = -^^-, r 3 3 = r 2 2 = ^cott 1 . (A.20) 

We can use the Maurer-Cartan forms to analyze the differential geometry of S 3 . The dreibein 
of S 3 is proportional to co a , and we have 

< = \ H • ( A - 21 ) 

In terms of forms, we have 

e a = e«dx» = r -u a . (A.22) 

Indeed, one can explicitly check that 

effiai, = G>. (A.23) 

The inverse vierbein is defined by 

E% = VabG^et (A.24) 
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which can be used to define left-invariant vector fields 



** = (A-25) 



Let us give their explicit expression in components: 



2 / d sin £3 d d 
£1 = - [ cost 3 — + -— - — - smt 3 cotti — 
r \ ot\ sin ti OT2 ots 

2 / . 9 costs 9 d . , - ... 

4 = - [smt 3 - — ^— +cost 3 cotti— , (A.26) 

r V ati sinti OT2 OT3 



4 

Of course, they obey 



r dts 



e a (4) = C (A.27) 
as well as the following commutation relations 

[4,4] = --€afc4- (A.28) 
r 

This can be checked by direct computation. If we now introduce the operators L a through 

4 = -La- (A.29) 
r 

we see that they satisfy the standard commutation relations of the SU{2) angular momentum 
operators: 

[L a ,L b ] = ie abc L c . (A.30) 
The spin connection uj\ is characterized by 

de a + uj a b A e b = 0. (A.31) 

Imposing no torsion one finds the explicit expression, 

«%, = (V, - rJ„eS) , (A.32) 

or, equivalently, 

d,e a v = T^el - e^V (A.33) 

In our case we find 

< = (A.34) 

A. 3 Laplace Beltrami operator and scalar spherical harmonics 

The scalar Laplacian on S 3 can be calculated in coordinates from the general formula 

or equivalently 

-A = Grdpdv - G^. (A.36) 
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In this case it reads 

- A ° = h (5 + cot *k + csc2 h m + csc2 tx M " 2 csc h cot 41 ^) ■ (A ' 37) 

It is easy to check that it can be written, in terms of left-invariant vector fields, as 

-A° = ^e (A.38) 

a 

To see this, we write 

E % = E E Xd,d v + E Kq£qJ- (A.39) 

a a ct 

The first term reproduces the first term in (A. 36). We now use the identity 

d,El = E^ - T^El (A.40) 

After contraction with E% and use of the explicit form of the spin connection, we see that only 
the second term survives, which is indeed the second term in (A. 36). 

The Peter-Weyl theorem says that any square-integrable function on S 3 ~ SU(2) can be 
written as a linear combination of 

Sf n , m,n = l,---,dj (A.41) 

where 

Sj : SU(2) ^ M djXd . (A.42) 

is the representation of spin j and dimension dj, and M djXd . are the inversible square matrices 
of rank dj. The function S™ n is just the (m,n)-th entry of the matrix. The eigenvalues of 
the Laplacian might be calculated immediately by noticing that, in terms of the SU(2) angular 
momentum operators, it reads 

A = 4jL 2 , (A.43) 
and since the possible eigenvalues of L 2 are 

+ j = 0,^,---, (A.44) 

we conclude that the eigenvalues of the Laplacian are of the form 

A, = ^i(i + 1), 3 = 0, (A.45) 

Notice that the dependence on r is the expected one from dimensional analysis. The degeneracy 
of these eigenvalues is 

d) = (2j + l) 2 (A.46) 
which is the dimension of the matrix M djXdj . 
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A. 4 Vector spherical harmonics 

The space of one-forms on S 3 can be decomposed in two different sets. One set is spanned by 
gradients of S™™, and it is proportional to 

Sf\T a )fuj a . (A.47) 

The other set is spanned by the so-called vector spherical harmonics, 

VjT, e = ±l, m = l,--- ,d j± i, n = !,■■■ ,d jT i, (A.48) 

see Appendix B of [8] for a useful summary of their properties. The e = ±1 corresponds to two 
linear combinations of the uj a which are independent from the one appearing in (A.47). The 
vector spherical harmonics are in the representation 

3±Lj*l) (A-49) 



2 1J 1 2, 

of SU{2) x SU{2). We will write them, as in [8], as V a , where 

a = (j,m,m' , e), (A. 50) 
and we will regard them as one-forms. They satisfy the properties 

Sv a = 0, *dV a = -e(2j + l)V a . (A.51) 

It follows that 

*d*dU Q = -A 1 U a = (2j + l) 2 U Q . (A.52) 

Their degeneracy is 

2d j+ id j _i=4j(2j + 2). (A.53) 

A. 5 Spinors 

Using the dreibein, we define the "locally inertial" gamma matrices as 

la = Efrp, (A.54) 

which satisfy the relations 

{la, lb} = 25 ab , [la, lb] = Zitabclc- (A.55) 

The standard definition of a covariant derivative acting on a spinor is 

V M = + \uflalb = ^ + lb]- (A.56) 

Using the commutation relations of the gamma matrices j a and the explicit expression for the 
spin connection (A. 34) we find 

=df, + ^-eabc^abd^ld = d^ + ^e^7 c 

' (A.57) 
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It follows that the Dirac operator is 



= -vfdp + l~ r = -h a E^ + A = -i 7 X + A (A.58) 



Let us now introduce the spin operators 



Sa = \la, (A.59) 



which satisfy the SU (2) algebra 

[S a , S b ] = ie abc S c . (Am) 
In terms of the S a and the SU(2) operators L a , the Dirac operator reads 

-Up = U*L ■ S + I V (A.61) 



r V 2, 

The calculation of the spectrum of this operator is as in standard Quantum Mechanics: we 
introduce the total angular momentum 

J = L + S, (A.62) 

so that 

4L • S = 2 (J 2 - L 2 - S 2 ) . (A.63) 

Since S corresponds to spin s = 1/2, and L to j, the possible eigenvalues of J are j ± 1/2, and 
we conclude that the eigenvalues of (A.61) are (we set r = 1) 

2 ((i ± iVi ± i + lWi + 1 )W 2 i +i , fo,+ , (A.64) 



-2j-\ for-, 



with degeneracies 



These can be written in a more compact form as 



A± = ±(n+-1, d±=n(n + l), n=l,2,--- (A.66) 
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